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Appendices 


This  appendix  section  contains  the  various  codes  that  were  used  to  do  the 
numerical  analysis  discussed  in  Section  I. 


Also  included  here  are  two  papers  which  were  prepared  to  conclude  the  work 
carried  out  during  an  earlier  period  of  the  Grant,  preceding  the  work  on  focused 
beams.  This  work  was  to  be  in  preparation  for  an  anticipated  later  effort  to 
merge  the  focused  beam  work  and  the  earlier  work  on  plate  excitation  by  bounded, 
but  non-focused  beams.  Unfortunately,  the  results  presented  here  could  not  be 
used  in  conjunction  with  the  results  presented  in  Section  I  due  to  the  decision 


APPENDIX  A 
LFOCUS.EXE 
USER’S  GUIDE 
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A.I.  INTRODUCTION 


This  appendix  outlines  the  steps  necessaiy  for  using  the  line  focus  software, 
1JFOCUS.EXE,  developed  to  support  ONR  research  on  (^tical  probing  of  ultrasonic  fields. 
Section  A.2  outlines  the  hardware  requirements  for  program  execution.  Since  the  software 
is  menu  driven  rather  than  batch  processed.  Section  A3  provides  a  short  user’s  guide  for 
the  menu  driven  interface.  Section  A.4  describes  the  output  file  structure  used  for  data  post 
processing. 

A2.  Hardware  Requirements 

The  minimum  hardware  configuration  required  to  operate  LFOCUS.EXE  consists 
of  a  computer  with  an  Intel  80386  Central  Processing  Unit  (CPU)  or  80486  CPU  with  the 
following  features: 

•  one  high  density,  SVz”  or  SVi"  floppy  disk  drive, 

•  PC  Disk  Operating  System  (DOS)  Version  3.3  or  later’, 

•  Monochrome  display  or  display  adapter, 

•  Intel  80387  Math  CO-Processor^, 

•  A  hard  drive  with  40  to  50  Mega-Bytes  (Mb)  of  available 
storage,  and 


Program  ccMnpatibility  has  been  verified  with  MS-DOS  v  3.3  and  4.01  and  Digital  Researdi 
DOS  V  5.0  and  6.0. 
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The  80486  chip  has  an  integral  co-processor. 


4.0  Mb  of  extended  memory. 


LFOCUS.EXE  accesses  the  Intel  co-processor  and  wiU  not  run  without  R.  As  an 
alternative,  the  program  can  be  compiled  to  access  the  Weitek,  ABACUS  chip.  This  change 
could  improve  program  execution  speed  by  about  a  factor  of  two.  A  graphics  monitor  does 
not  improve  program  performance.  Additional  hard  disk  storage  ^ace  is  only  useful  for 
multiple  program  runs.  A  Microsoft®  or  compatible  mouse  can  be  used  as  a  prating 
device  for  program  execution;  however,  a  mouse  is  not  required  for  program  execution. 

A3.  Initial  Execution 

To  run  LFOCUS.EXE  for  the  first  time  change  to  the  hard  disk  directory  which 
contains  the  program,  type  LFOCUS,  and  press  <Enter>.  LFOCUS.EXE  loads  and 
presents  the  initial  menu  screen  shown  in  Figure  A.l. 


Data  Entry  Window 


Fib  Opens  Menu  <E8C>  Closes  Menu  Alt-X  Exits  Prograa 


Figure  A.l.  LFOCUS.EXE  Initial  Menu 


If  a  mouse  is  present,  move  the  mouse  cursor  over  the  desired  functicm  and  click  the 
left  mouse  button  to  execute  the  function.  If  a  mouse  is  not  present,  press  the  highlighted 
key  to  execute  the  function.  For  example,  the  initial  menu  can  be  activated  1^  pressing  the 
Alt  key  and  W  to  select  the  Data  Entry  Window,  moving  the  mouse  cursor  over  Data  Entry 
Hindow  and  clicking  the  left  mouse  button,  or  pressing  FIO.  If  the  Data  Entry  Window 
functi(Mi  is  executed  by  either  teclmique,  the  first  pull-down  menu  appears  as  shown  in 
Figure  K2. 


Figure  A.2.  Data  Entry  Options 


Run  Parameters  (F6)  specifies  the  number  of  spectra  used  in  the  calculations 
(minimum  and  maximum)  and  the  on  and  off-axis  data  output  points.  Off-Axis  Data  Points 
(F7)  specifies  the  on-axis  locations  wdiere  the  off-axis  data  points  are  output  to  a  disk  file. 
Transducer  (F8)  specifies  the  key  parameters  associated  with  the  line  focus  transducer,  such 
as,  dimensions,  aperture  angle,  initial  on-axis  pressure,  etc..  Boundary  Values  (F9)  provides 
options  for  specifying  the  pressure  amplitude  distribution  at  the  transducer  surface.  To 
activate  a  menu  press  the  highlighted  key,  the  specified  function  key,  or  use  the  mouse 
cursor  and  click  on  the  choice. 


A.4 


A.3.1.  Run  Parameters 


To  specify  the  input  parameters  for  a  run  press  R,  function  key  F6,  or  highlight  the 
Run  Parameters  option  with  the  mouse  cursor  and  clicL  The  Run  Parametm  menu  shown 
in  Figure  A3  appears.  The  parameter  values  displayed  are  those  used  the  last  time  the 
program  was  executed.  To  exit  this  menu  without  changing  any  values,  use  the  mouse 
cursor  to  click  the  [■]  symbol  or  press  <Enter>.  To  move  through  the  parameters  on  the 
menu  without  a  mouse,  press  the  Tab  key.  Each  time  the  Tab  key  is  pressed  the  next  menu 
option  is  highlighted  for  data  entry.  To  back-up,  press  Shift  Tab  and  the  previous  data 
entry  option  is  highlighted  for  data  entry.  To  select  a  field  for  data  entry  move  the  mouse 
cursor  over  the  text  describing  the  value  and  click. 


Data  Entry  Window 


Run  Paraaeters 


Uaa  Tab  Kay  or  UMiaa  to  aalact  data  antry  flald 

UlnlBua  #  Spectra 
Maxiaua  »  Spectra 
Nuaber  On-Axis  Points 
Nuaber  Off-Axis  Points 
Radial  Step  Size 
Min  On -Axis  Output 
Max  On-Axis  Output 
Max  Off-Axis  Range 
Max  Off-Axis  Points 

Error  Message 


(4  to  49): 

5 

(Min  »  to  50): 

26 

(100  to  400): 

300 

(1  to  40): 

30 

(0.02  to  0.001): 

0.005 

(0.001  to  0.99); 

0.001 

(1.05  to  8.0): 

1 .4 

(1.1  to  10.0): 

10 

(10  to  200): 

200 

F10  Opens  Menu  <Esc>  Closes  Menu  Adt-X  Exits  Progri 


Figure  A3.  Run  Parameters  Menu 


To  enter  in  a  new  value  key  in  the  appropriate  number.  Do  not  press  Enter,  if  you 
wish  to  change  other  values  before  exiting  the  menu.  Pressing  <Enter>  will  exit  the  menu, 
if  there  are  no  errors.  If  a  value  is  keyed  in  which  is  out  of  range,  an  enor  message  appears 
when  you  attempt  to  exit  the  menu.  For  example.  Minimum  #  Spectra  represents  the 
number  of  spectral  terms  present  at  the  start  of  a  calculation.  As  shown,  this  is  allowed  to 


take  on  a  value  between  4  and  49  terms.  If  a  value  CHitside  of  this  range  is  keyed  in,  such 
as  2,  the  program  writes  the  error  message  shown  in  Figure  A.4  when  the  <Enter>  key  is 
pressed.  The  program  will  not  exit  the  menu  until  the  entry  is  corrected.  If  multiple  errors 
are  made  when  keying  in  values,  each  enor  has  to  be  corrected. 


Data  Entry  WiiMlow 


flun  Paraaeters 


Um  Tab  Kay  or  Mouaa  to  aalact  data  antry  flald 

mnlwa  t  Spaotra 
Maxiaua  #  Spectra 
Nuaber  On-Axis  Points 
Nuabar  Off-Axis  Points 
Radial  Step  Size 
Uln  On -Axis  Output 
Uax  On-Axis  Output 
Max  Off-Axis  Range 
Uax  Off-Axis  Points 

Error  Message  lUniaua  Spectral  Taro  Out  Of  Rangai 

Ok  ■  Cancel 


(4  to  4»): 

2 

(Min  f  to  50): 
(100  to  400): 

26 

300 

(1  to  40): 
0.02  to  0.001): 

30 

0.005 

0.001  to  0.99): 

0.001 

(1.05  to  2.0): 

1.4 

(1.1  to  10.0): 

10 

(10  to  200): 

200 

PiO  Opens  Menu  <Esc>  Closes  Menu  Alt'X  Exits  Prograa 


Figure  A.4.  Typical  Error  Message 
The  data  entry  values  covered  in  the  Run  Parameters  menu  include: 


Minimum  #  Spectra 
Maximum  #  Spectra 

Number  On-Axis  Points 

Number  Off-Axis  Points 

Radial  Step  Size 


minimum  number  of  spectral  terms  present  at  the 
start  of  a  calculation, 

maximum  number  of  spectral  terms  allowed  in 
the  calculation,  that  is  the  number  of  terms 
allowed  before  the  Fourier  spectra  is  truncated, 

number  of  on-axis  points  where  the  amplitude 
and  phase  spectra  are  saved  during  the 
calculation, 

number  of  points  on  the  axis  i^ere  the  off-axis 
amplitude  and  phase  spectra  are  saved  during  the 
calculation, 

normalized,  transverse  mesh  step  size  used  in  the 
calculation. 
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the  oD-axis  point  where  the  oo-axis  output  starts 
(0  is  at  the  face  of  the  transducer  and  1  is  at  the 
geometric  focal  point), 

the  on-axis  pcant  where  the  calculation  is 
terminated  and  the  last  on-axis  data  values  are 
written  to  a  file, 

maximum  normalized  size  the  transverse  grid  is 
allowed  to  expand  to  during  the  calculations,  and 

maximum  number  of  off-axis  points  which  are 
output  at  each  ^cified  output  location. 

The  values  shown  in  Figures  A3  and  A.4  within  the  brackets  are  the  allowable  data  enby 
ranges. 

A3.2.  Off-Axis  Data  Point  Menus 

To  review  or  change  the  on-axis  locations  selected  for  off-axis  data  output  prior  to 
running  LFOCUS.EXE,  select  the  data  entiy  window  and  press  O,  F7,  or  use  the  mouse 
cursor  to  select  Off-Axis  Data  Points  from  the  cations  shown  in  Figure  A.2.  At  this  point, 
the  Off  Axis  Location  1  menu  shown  in  Figure  A3  appears.  Use  the  Tab  or  Shift  Tab  keys 
to  move  between  the  data  entry  fields.  To  exit  without  changing  any  values  use  the  mouse 
to  click  on  the  [■]  symbol  or  the  Ok  or  Cancel  buttons.  Alternatively,  use  the  Tab  or  Shift 
Tab  keys  to  move  to  the  Ok  or  Cancel  buttons.  The  range  of  allowable  off-axis  data  values 
is  determined  1^  the  Run  Parameters  menu.  For  examjrfe,  in  menu  shown  in  Figure  A3, 
the  maximum  on-axis  output  point  is  1.4  (40%  beyond  the  geometric  focal  point  for  the 
transducer).  Data  entry  follows  the  techm'que  previously  described  for  the  Run  Parameto^ 
menu.  Remember,  the  software  will  not  exit  a  menu  which  has  a  value  out  of  range.  Also, 
the  program  expects  the  values  to  be  monotonicaOy  increasing.  This  menu  will  continue  to 
appear  until  all  points  have  been  selected  for  off-axis  data  output.  Once  the  last  data  point 
has  been  selected,  the  program  returns  to  the  initial  data  entry  menu  shown  in  Figure  A.I. 


Min  On-Axis  Output 

Max  On-Axis  Output 

Max  Off-Axis  Range 
Max  Off-Axis  Points 
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P10  opens  Menu  <Eso  Closes  Menu  Alt-X  Exits  Prograa 

Figure  AS.  Off-Axis  Data  Point  Selection 
A3.3.  Transducer  Qiaracteristics  Menu  Data  Entry 

To  review  or  change  the  transducer  data  values  prior  to  running  U'OCUS.EXE  select 
the  data  entry  window  and  press  T,  F8,  or  use  the  mouse  cursor  to  select  IVansducer  from 
the  options  shown  in  Figure  A.2.  At  this  point  the  Transducer  Characteristics  menu  shown 
in  Figure  A.6  appears.  Use  the  Tab  or  Shift  Tab  keys  to  move  between  the  data  entry 
fields.  To  exit  without  changing  any  values  use  the  mouse  to  cb'ck  on  the  [■]  symbol  or  the 
Ok  or  Cancel  buttons.  Alternatively  use  the  Tab  or  Shift  Tab  keys  to  move  to  the  Ok  or 
Cancel  buttons.  Data  entry  follows  the  technique  previously  descnhed  for  the  Run 
Parametors  menu.  Remember,  the  software  will  not  allow  the  user  to  exit  a  menu  ^hich 
has  a  value  out  of  range.  Once  the  last  data  entry  has  been  made  the  program  returns  to 
the  initial  data  entry  menu  shown  in  Figure  A.I. 


Data  Entry  Wirwlo* 


Um  Tab  My  or  Hwmo  to  ooloct  tfoto  ontry  floM 


Oiaaeter  (Haters:  0  ..  1)  0.023009 

rc.al  Length  (Haters:  0  ..  1)  0.04445 

Initial  Frequency  (mz:  0  ..  10)  1.0M7 

Initial  Pressure  (kPa:  0  ..  2000)  57.028 


Figure  A.6.  Transducer  Data  Entry 


AJ.4.  Boundary  Value  Menus 


A  three  level  menu  system  is  used  for  selectii^  an  existing  or  entering  a  new 
normalized  pressure  distribution  across  the  face  of  the  Une  focus  transducer.  To  select  this 
option  from  the  choices  presented  in  Figure  A.2,  press  &  F9,  or  highlight  Boundary  Value 
with  the  mouse  cursor  and  click.  The  Transducer  Pressure  Distribution  Options  Menu 
shown  in  Figure  A.7  appears.  This  menu  provides  three  choices  for  the  pressure 
distribution  at  the  transducer  surface. 


Previous;  The  distribution  used  in  the  previous  run 
Uniform;  A  uniform  distribution,  e.g.  the  peak  pressure  at 
each  location  on  the  surface  of  the  transducer  isa 
assumed  to  be  the  same 

New;  The  pressure  distribution  will  be  entered  from 

the  key  board  or  a  data  file 
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Data  Entry  Windon 


[■ITransducer  Pressure  Distribution  Options 

To  Select  Option,  Press  HighliQbted  Letter 

Pressure  Distribution  Source 

(•)  Previous  (•)  Use  Pile 

(  )  Unlfora  (  )  Enter  Data 


P10  Opens  Menu  <Esc>  Closes  Menu  Alt*X  Exits  Prograa 


Figure  A.7.  Transduce  Pressure  Distributi(»  Options 


To  select  a  vahie,  either  press  the  highlighted  letter  or  move  the  mouse  cursor  over  the 
option  and  click.  The  default  is  to  use  the  normalized  pressure  distribution  for  the  previous 
run.  To  change  to  uniform  distribution,  press  U  or  mo^e  the  mouse  cursor  over  Uniform 
and  click.  To  enter  a  new  pressure  distribution,  press  N  or  use  the  mouse  cursor  to  click 
on  New.  If  New  is  selected,  the  source  choices  become  active.  The  default  is  to  read  an 
existing  data  file  which  describes  the  pressure  distributicMi.  To  enter  a  new  data  file,  press 
E  or  click  on  Enter  Data  with  the  mouse  cursor.  Once  the  final  choice  has  been  made, 
press  <Enter>,  O,  or  click  on  Ok  with  the  mouse  cursor. 

If  the  source  option  is  Use  File,  the  Pressure  Distribution  File  Name  Menu  shown 
in  Figure  A.8  appears.  Press  F,  use  the  Tab  key,  or  the  mouse  cursor  to  select  File  Name. 
Key  in  the  file  name  and  press  <Enter>.  If  an  error  occurs  during  the  data  entiy,  a 
descriptive  message  will  appear  in  the  Error  Message  window. 
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Figure  A-8.  Entering  New  Ffle  Name 


Once  the  desired  field  name  has  been  selected  press  <Enter>.  To  terminate  the  process 
at  any  time,  tab  to  the  Cancel  button  or  click  on  this  button  with  the  mouse  cursor. 

A  simQar  process  is  followed  in  generating  a  new  data  file.  Select  the  Enter  Data 
option  from  the  menu  shown  in  Figure  A.7  and  press  <Enter>.  The  New  Pressure 
Distribution  File  Moiu  shown  in  Figure  A.9  appears.  Press  F,  use  the  Tab  key,  or  the 
mouse  cursor  to  select  File  Name.  Key  in  the  file  name.  Use  the  Tab  key  or  the  mouse 
cursor  to  select  Number  of  points.  Key  in  a  value  between  1  and  100  and  press  <Enter>. 
If  an  error  occurs  during  data  entry,  a  descriptive  message  appears  in  the  Error  Message 
window.  Use  the  Tab  key  or  the  mouse  cursor  to  move  to  the  data  entry  window  that 
contains  the  error  and  key  in  a  corrected  value.  To  exit  at  any  time  without  entering  a  new 
data  file,  use  the  Tab  key  or  the  mouse  cursor  to  select  the  Cancel  butt  n  and  press  enter 
or  click  on  the  button. 


Data  Entry  Window 


New  Pressure  Distribution  File  Naae 


Pile  Naae  (<40  Char) 
Nuaber  of  points 


P10  Opens  Menu  <Eso  Closes  Menu  Alt>X  Exits  Prograa 


Figure  A.9.  Entering  a  new  file  name  and  number  of  points 


Data  Entry  Window 


—  ■  Pressure  Distribution  Location  -  1  — 

To  Change  The  Value,  Press  Highlighted  Letter 

Radial  Point  (a)  0.0 

Noraalized  Pressure  1.0 


Error  Message 


Cancel 


Figure  A.10.  Altering  New  Pressure  Distributitn  Ffle 

To  enter  the  transducer  amplitude  distribution  use  the  menu  shown  in  Figure  A.10. 
First  use  the  Tab  key  or  the  mouse  to  select  Radial  Point  (m).  Ihe  radial  location  is 
measured  normal  to  the  transducer  axis  of  symmetry.  The  first  entry  must  be  the  value  on 
the  ^rmmetry  axis,  e.g.  0.  After  keying  in  the  appropriate  value  use  the  Tab  key  or  the 
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mouse  to  select  Normalized  Pressure.  Key  in  a  value  between  0  and  1.  If  the  radial 
location  extends  beyond  the  transducer  boundary,  an  error  message  appears.  Also,  a 
normalized  pressure  greater  than  unity  generates  an  error  message.  After  the  location  and 
pressure  value  have  been  correctly  entered,  press  <Enter>  or  the  Ok  button.  After  the  last 
data  entry  is  made,  the  program  returns  to  the  program  main  menu  shown  in  Figure  A.1. 

A3.5.  Program  Execution 

After  data  entry  is  completed,  press  Alt-X  or  use  the  mouse  cursor  to  click  the  Alt-X 

Exits  Program  section  of  the  menu  shown  in  Figure  A.1  to  start  program  execution. 

Depending  on  the  data  values  used  and  computer  CPU  speed,  program  execution  can  run 

from  several  hours  to  as  much  as  a  couple  of  days.  While  the  program  is  executing,  the 

status  message  shown  below  is  written  to  the  screen. 

Nstep  =  264  z  =  0.136237  k  =  29  m  =  4  Imax  =  113 

Nstep  =  273  z  =  0.140900  k  =  30  m  =  4  Imax  =  115 

Nstep  =  282  z  =  0.145563  k  =  31  m  =  4  Imax  =  115 

Nstep  is  the  calculation  step  number.  Z  is  the  normalized  on-axis  distance  from  the  face 

of  the  transducer.  K  is  the  number  of  the  output  point.  The  current  number  of  harmonic 

terms  being  used  in  the  calculation  is  specified  by  m.  Finally,  Imax  is  the  current  location 

of  floating  outer  boundary  for  the  edge  of  the  calculational  mesh. 

To  interrupt  program  execution  press  Ctrl-C  or  Ctrl-Break.  After  two  screen  writing 
^cles,  all  relevant  data  is  saved  to  the  hard  disk  and  normal  program  termination  occurs. 
To  restart  the  program,  type  LFOCUS  and  press  <Enter>.  When  the  program  has  finished 
running,  it  writes  PROGRAM  EXECUTION  COMPLETED  to  the  screen. 

A.4.  Output  Data  File  Structures 

After  LFOCUS.EXE  has  completed  execution  four  text  data  files  are  written  - 
amplitude  and  phase  files  for  on-axis  and  off-axis  data  points.  Figure  A.ll  summarizes  the 
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relevant  characteristics  of  AMPZ  and  PHASEZ  -  the  on-axis  amplitude  and  phase  spectra 
output  files.  The  first  44  lines  in  each  file  summarizes  the  input  data  used  to  conduct  the 
run.  This  is  followed  in  each  case  by  the  amplitude  spectra  or  phase  spectra,  depending 
upon  the  file  type.  Since  the  number  of  spectral  terms  will  vary,  depending  upon  the 
nonlinearity  of  problem,  e.g.  initial  pressure,  effective  focusing  gain,  etc.,  the  first  output  is 
the  number  of  ^ctral  terms,  the  second  output  is  the  location  along  the  z  axis.  The 
subsequent  output  terms  are  either  the  amplitude  or  phase  spectra  depending  upon  the  file 
type.  The  amplitude  ^ctral  terms  are  normalized  to  the  peak  pressure.  The  phase  terms 
are  in  radians. 


Aapz.  Data  Structure 
Lines  1  through  44  >  file  header 

Line  45  through  end  «  aaplltude  spectra  with 
■  2  Aid  Ain 

4  S.6633E-0003  O.OOOOE^OOOO  1.0032E-HMOO 

4  1.0327E-0002  O.OOOOE+OOOO  1 .0057E+0000 

4  1.4990E-0002  O.OOOOE-KXKK)  1  .OOSIE'HXKX) 


AI21  ••• 
6.2014E-0005  ••• 
1.0892E-0004  ••• 
1.5616E-0004  ••• 


Phasez.  Data  Structure 
Lines  1  through  44  >  file  header 

Line  45  through  end  »  phase  spectra  with 

■  z  PhlfOl  Phil  11 

4  5.6633E-0003  O.OOOOE^OOOO  •4.8870E-0006 

4  1.0327E-0002  O.OOOOE+OOOO  -1 .460SE-0005 

4  1.4990E-0002  0.0000£*0000  •2.96B9EOOOS 


Phl|2| 

•6.4645E  0006 
-1 .8326E-0005 
-3.6417E-0005 


•  •• 
•  •• 
•  •• 
•  •• 


Figure  A.ll.  Ampz  and  Phasez  Output  Data  Structure 


A  simQar  file  structure  is  used  for  the  off-axis  points,  as  shown  in  Figure  A.  12. 
Again,  the  first  44  lines  of  AMPZR  or  PHASEZR  summarize  the  input  data.  Line  nine  is 
an  integer  entry  wfaidi  ^cifies  the  number  of  off-axis  data  output  points.  Starting  with  the 
first  output  location,  a  record  is  written  which  contains  the  number  of  spectral  terms  used 
at  the  location,  the  normalized  on-axis  location,  the  normalized  off  axis  location,  and  the 
first  through  n*^  ampfa’tude  or  phase  terms. 


A. 14 


Lines  1  through  44  ^ 
Line  9  * 

Line  45  through  end  = 

■  z 

5  7.5000E-0001 

5  7.SOOOE0O01 
5  7.SOOOE-0OO1 


Aapzr.  Data  structure 
file  header 

nuaber  of -axis  data  output  locations 
aaplitude  spectra  with 

r|ll  ACOl  AMI 

O.OOOOE^OOOO  O.OOOOE+OOOO  2.833^^0000 

rI2I 

2.0000E-0002  O.OOOOE^OOOO  2.7830E>0000 


rl31 

4.0000E-0002 


0.0000540000  2.8362E40000 


•  •• 

•  •• 

•  •• 


Lines  1  through  44  • 
Line  9  * 

Line  45  through  end  * 

■  z 

5  7.5000E-0001 

5  7.5000E-0001 

5  7.5000E-0001 


Phasezr.  Data  Structure 
file  header 

nuaber  of-axls  data  output  locations 
aaplitude  spectra  with 


r|l) 
0.0000E40000 
rl2I 

2.0000E-0002 


4.00^- 


0002 


Phl(Oj 
O.OOOOE40000 

O.0000E4O00O 

O.000OE4O00O 


9.7M7IU0OO 

1 .0189E40001 
1 .1529E40001 


•  •• 
•  •• 

•  •• 

•  •• 


Figure  A.12.  Ampzr  and  Phasezr  Output  Data  Structure 


The  amplitude  spectra  are  normalized  relative  to  the  initial  input  pressure  and  the  phase 
spectral  terms  are  in  degrees. 
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APPENDIX  B 

MENU.PAS 
PROGRAM  LISTING 


© 
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This  appendix  provides  a  tisting  of  the  menu  program  used  by  LFOCUS.EXE.  This 
program  is  written  in  Turbo  Pascal  v  6.0  and  called  by  the  FORTRAN  code  during  initial 
execution.  Since  the  program  uses  the  standard  Turbo  Vision  applications  framework,  only 
the  unique  source  code  is  commented. 


prograa  MENU; 

use«  Obiects, 
Drivers , 
Views, 
Menus, 
Dialogs, 
App; 


const 

MaxLlnes  >  100; 

WlnCount:  Integer  >  0; 

osFlleOpen  ■  100; 

caNesWln  ■  101; 

cMewOlalog  >  102; 

caOffAxlsOata  >  103; 

caNewTransducer  *  104; 

caRunInput  *  105; 

caBoundary  »  106; 

caPrevloua  >107; 

caUnlfora  *  108; 

caNeaFlle  ■  109; 

caNewFlleValues  *  110; 

var 

LlneCount , 

MaxNuaber, 

OffAxlsIndex, 

code, 

n  :  INTEGER; 

Lines:  array(0. .MaxLlnes  -  1]  of  PStrlng; 
errorflag  :  BOOLEAN; 


type 

{  Record  types  used  by  ibe  program  for  temporary  data  I/O 


RunlnputValues 

MINM 

MAXH 

Npolnts 

NOpolnts 

OeltaR 

Zaln 

Zaax 

Rgaax 

RadlusMax 

erroraessage 

end; 


record 


:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

:  STRING 

128 

TransducerOata 

Dlaaeter 

RadlusOfC 

Frequency 

Pressure 

erroraessage 

end; 


RECORD 

STR1N6M28 

string} 128 
string} 128 

8TRIN6I128 

string} 128 


ZuOata  >  RECORD 
location  :  STR1NGI128I; 
erroraessage  :  STRING} 128]; 
end; 


BoundaryData  >  RECORD 
RadloButtonDataA  :  Word; 
RadloButtonOataB  :  Word; 
end; 


FlleNaaeData  ^  RECORD 
Naae  :  STRING} 128]; 

ErrorMessage  :  STRING} 128]; 
end; 
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FlleNaaeOataa  «  RECOno 
Naae  :  STRINGI128]; 

NuaberPoints  :  STRING{l28i; 
Errorilessage  :  STRIMG(128i; 
end; 

BoundaryOataValues  -  RECORD 
location  :  STRINGM28]; 
NoraalizedP  :  STRING ( 1281: 
erroraessage  :  STRING(l28l: 
end; 


TVyApp  >  ob)ect(TAppllcatlon) 

procedure  HandleEvent(var  Event:  TEvent);  virtual; 
procedure  InitMenuBar:  virtual; 
procedure  InitStatuaLine;  virtual; 
procedure  Nea^ransducer; 
procedure  RunParaaetera; 
procedure  OffAxlaOata; 
procedure  NeeMlndow; 
procedure  SelectBoundary; 
procedure  EnterFlleNaae; 
procedure  EntarFlleNaae2 : 
procedure  BoundaryAxlsOata; 
end; 

Plnterlor  »  ‘Tlnterlor; 

Tlnterlor  =  object (TScroller) 

constructor  Inlt(var  Bounds:  TRect;  AHScrollBar,  AVScrollBar:  PScrollBar); 
procedure  Draw;  virtual; 
end; 

POeaoMrindow  «  ''TOeaoWlndow; 

TDeaoWlndow  =  object (TWlndow) 

Rlnterlor,  Llnterlor:  Plnterlor; 

constructor  Ir.it(Bounds;  TRect;  WlnTitle:  String;  WindoeiNo:  Word); 
function  Uakeinterlor (Bounds:  TRect;  Left:  Boolean):  Plnterlor; 
procedure  SlzeLlalts(var  Min,  Max:  TPolnt);  virtual; 

end; 

POeaoOlalog  *  ‘'TOeaoOlalog; 

TDeaoOlalog  -  object (TOlalog) 
end; 


var 

Transducer 

Runinput 

Zu 

BoundaryCholce 

FlleNaae 

FlleNaae2 

BoundaryCholceValues 

NuaberSoundaryPolnts 

textflle 


TransoucerOata; 

Runinput Values: 

Array(0..50]  of  ZuOATA; 

BoundaryOata; 

FlleNaaeOata; 

FlleNaaeOataZ; 

ArraylO. .100]  of  BoundaryOataValues; 
INTEGER; 

TEXT; 


constructor  Tlnterlor. I nit (var  Bounds:  TRect;  AHScrollBar,  AVScrollBar:  PScrollBar); 
begin 

TScroller. Inlt (Bounds,  AHScrollBar,  AVScrollBar); 

Options  :=  Options  or  ofFraaed; 

SetLlBlt(l28,  LlneCount); 
end; 
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procedure  Tlnterlor.DroM; 
var 

Color:  Byte; 

1,  Y:  Integer; 

B:  TOrawBuffer; 
begin 


eoin 

Color  GetColor(l); 

for  Y  :»  0  to  Slze.Y  -  1  do 

begin 

fioveChar(B,  '  Color,  Size.X); 


i  :»  Delta. Y  ♦  Y: 
if  (I  <  LineCount)  and  (Lineal I ]  <>  nil)  then 
lloveStr(B,  Copy(Line8(I|‘,  Delta. X  *  1,  Size.X),  Color); 

WriteLine(0,  Y,  Size.X,  1,  B); 
end; 
end; 

conatructor  TDeaoWlndo«.lnit(Bounda:  TRect;  WinTitle:  String;  wlndcwNo:  ttord) 
var 

S:  atringIS]; 

R:  TRect; 
begin 

8tr(WindoiaNo,  8); 

TWindow.lnitiBounda,  WlnTltle  *  '  ‘  *  S,  wnNoNuaber); 

GetExtent (Bounda) ; 

R.A8aign(Bounda.A.X,  Bounda. A. Y,  Bounda. B.X  div  2*1,  Bounda. B.Y); 
LInterior  :■  llakeInterior(R,  True); 

Llnterior'.GroeMode  gfQroMHiY; 

I naertf LInterior) ; 

R.Aaalgnl Bounda. B.X  div  2,  Bounda. A. Y,  Bounda. B.X,  Bounda. B.Y); 

RInterlor  :=  Makelnterlor(R,Fal8e) ; 

Rlnterior'.GroaMode  gfGrowHlX  *  gfGrcHdliY; 

Inaert(RInterior); 

end; 

function  TDeaoWindow.Makelnterlor (Bounds:  TRect;  Left:  Boolean):  Plnterior; 
var 

HScrollBar,  VScroIlBar:  PScrollBar; 

R:  TRect; 
begin 

R.A88ign(Bound8.B.X-1 ,  Bounds. A. Y-f  1 ,  Bounds. B.X,  Bounds. B. Y-1 ); 

VScrollBar  :*  Hem  (PScrollBar,  Init(fl)); 

VScroIlBar* .Options  :=  VScrollBar* .Options  or  ofPostProcess; 
if  Left  then  V8croll8ar*.Gro«Mode  gfGrowHiY; 

InaertlVScrollBar) ; 

R.Asslgn(Bounds.A.Xe2,  Bounds. B. Y-1 ,  Bounds. B.X-2,  Bounds. B.Y); 

HScrollBar  :»  New(PScrollBar .  Init(R)): 

HScrollBar* .Options  HScrollBar*. Options  or  ofPostProcess; 
if  Left  then  HScrollBar *.Gro«illode  :=  gfGroeHiY  gfGroeLoY; 

Insert (HScrollBar ) ; 

Bounds. Grow(-1 ,-1 ); 

Hakelnterlor  :>  New(Plnterlor,  Init(Bounds,  HScrollBar,  VScrollBar)); 
end; 

procedure  TOeeoWlndow.SlzeLliits(var  Min,  Max:  TPolnt); 

var  R:  TRect; 

begin 

TWindow.SizeLlBlt8(Mln,  Max); 

Min.X  :*  LInterior* .Size.X  ♦  9; 
end; 
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This  is  the  event  handler  used  to  caD  the  data  entry  procedures. 


anNewTraosduoer 

omRunlnput 

onOSAxisDala 

cmNewWin 

onBoundaiy 


oaB  to  NeWTransducer 
call  to  RunParamctcrs 
can  to  OQAjdsDala 
can  to  NewWiodow 
can  to  SelectBoundary 


procedure  T1lyApp.HandleEvent(var  Event:  TEvent); 
begin 

TAppllcatlon .HandleEvent (Event ) : 

If  Event .Miat  ■  evCoaeana  then 
begin 

case  Event. Coaaand  of 

caNeuMfln  :  NewWindow; 

caNewTransducer  :  Neifrransducer; 
caftunlnput  :  RunParaaeters; 

caOffAxisData  :  begin 
VAL(Runlnput  .MDpointa  iMaxNuaber  .Code ) ; 

For  OffAxislndex  :•  1  To  MaxNtiaber  Do  OffAxisData: 
end; 

caBoundary  :  begin 

SelectBoundary: 

Case  Boundarycnoice.RadioButtonOataA  Of 
0  : ; 

1  :  begin 

{This  section  sets  a  unifora  aperture  distribution} 

NuaberBoundaryPoints  l ; 

BoundaryChoiceValues(0| .location  :=  ’O’: 

BoundaryChoiceValuesjoj .NoraalizedP  :s  ’1.6’; 
end; 

2  :  Case  BoundaryChoice.RadloButtonOataB  OF 

0  :  begin 

EnterFileNase; 

IF  skipflag  »  FALSE  THEN 
begin 

Assign(textf ile.FileNaae.naae) ; 

RESET ( textf ile); 

REAOLN ( textf ile .NuaberBoundaryPoints ) ; 

FOR  n  :=  0  TO  NuaberBoundaryPoints  •  1  00 
begin 

REAOLN(textfile,BoundaryChoiceValues[n)  .location) ; 
REAOLN(textrile,BoundaryChoiceValuesinj .NoraalizedP) ; 
BoundaryChoiceValues[n]  .erroraessage  :-  ”; 
end; 

FOR  n  :-  NuaberBoundaryPoints  TO  100  DO 
WITH  BoundaryChoiceVaiues(nj  DO 
begin 
location 
NoraalizedP 
erroraessage  :=  ” ; 
end; 

CLOSE(textfile); 

end; 

end: 

1  :  begin 

EnterFileNaae2: 

IF  skipflag  >  false  THEN 
begin 

VAL  ( FileNane2 .  NuaberPoints  .NuaberBoundaryPoints  .Code ) ; 

For  OffAxislndex  :«  0  to  NuaberBoundaryPoints  •  1  Do  BoundaryAxisOata; 
end; 
end; 

end; 

end; 

end 

else 

Exit; 

end; 

ClearEvent ( Event ) ; 
end; 
end; 
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This  is  the  interface  for  the  puU-down  menu  system,  ft  connects  the  commands  »ith  the  keystrokes  that  activate 
each  individual  command.  For  eiample.  R  or  Ft  are  linked  to  cmRun  Input  which  is  used  to  call  the  piooeihire 


procedure  TUyApp.InltHenuBar; 

var  R:  TRect; 

begin 


begin 

GetExtentfR); 

R.B.Y  :»  R.A.Y  ♦  1; 


HenuBar  :=  NewfPMenuBar,  Init(R,  NesMenut 

NewBubHenuf ’Data  Entry  ~W~lrKkw‘ ,  hcNoContext,  NewMertuf 
Ne«Itee('~R~un  Paraaetera* , 'F6*  ,kbF6,  caRunInput,  hcMoContext, 

Newlteai 'Axis  Data  Points’ ,'F7‘ ,kbF7,  caOffAxisData,  hcNoContext, 
Ne«ltea(‘~T~ransducer’ .  ‘F8*,  kbrB,  caNestTransducer,  hcNoContext, 

NewItea(’~B*oundary  Values’ ,  ‘F9’,  kbFg,  caBoundary,  hcNoContext, 

nil))))). 

nil)) 

end; 

{  This  prooeihire  controls  the  status  line  at  the  bottom  of  the  main  menu.  It  links  the  keys  that  open  and  dose  the 

main  menu. 

procedure  TMyApp.InitStatusLine; 

var  R:  TRect; 

begin 

GetExtent (R) ; 

R.A.Y  :*  R.B.Y  •  1; 

StatusLine  :>  NewfKtatusLine,  Inlt(R, 

NewStatua0ef(0,  SFFFF, 

NewStatutKey(‘*FlO~  Opens  Menu  <Esc>  Closes  Menu’,  kbFiO,  csHenu, 

NewStatusKey( ’"Alt-X”  Exits  Prograa’,  kbAltX,  caOuit, 
nil)), 
nil) 

V' 

end; 

{  Thu  procedure  is  used  to  check  real  input  values  to  determine  whether  they  are  valid.  If  the  dau  entry  is  not  a  real 

number,  the  error  flag  is  set  to  FALSE,  the  error  mesage  is  set  to  the  string  enorl.  and  the  procedure  branches  to 
the  exit.  If  the  dau  type  is  correct,  the  procedure  checks  to  see  whether  the  dau  entry  is  between  the  upper 
(Upperlimit)  and  lower  (LowerLimit)  limits.  If  it  is  not.  the  error  flag  is  set  to  FALSE,  the  error  message  is  set 
to  t^  string  enor2,  and  the  procedure  branches  to  the  exit 

procedure  CheckRealOatafVAR  NewVAlue,  NewValueZ  :  STRING; 

errorl ,  error2  :STR1HG; 

VAR  errorflag  :  BOOLEAN; 

LowerLinit,  UpperLinit  :  SINGLE); 

LABEL  localout; 

VAR  real value  :  SINGLE; 


Val(  NewVAlue,  real  value,  code) ; 

NewValueE 
IF  code  o  0  THEN 
begin 

NewValue2  :*  error! ; 

Errorflag  :«  FALSE; 

GOTO  localout; 
end; 

IF  (real value  <  lowerliait)  OR  (real value  >  upperlieit)  THEN 
begin 

Errorflag  :*  FALSE: 

NewValtie2  :»  errorz; 

GOTO  localout; 
end; 

localout  :; 
end; 
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{  This  piDoedure  b  used  to  check  toieger  input  values  to  determiae  whetber  they  are  valid.  If  the  data  entry  is  not 

an  integer,  the  error  flag  is  set  to  FALSE,  the  error  mesage  is  set  to  the  string  errorl,  and  the  procedure  branches 
to  the  ent.  If  the  data  type  is  correct,  the  procedure  cb^s  to  see  whether  the  data  entry  is  between  the  upper 
(UpperLonit)  and  lower  (LowerLimii)  limits.  If  it  is  not,  the  error  flag  is  set  to  FALSE,  the  error  message  is  set 
to  tte  string  ertor2,  and  the  piocerkire  branches  to  the  exii.  ) 

procedure  ChecklntegerOatafVAR  NewVAlue,  NewValue2  :  STRING; 

errorl,  error2  :STR1NG; 

VAR  errorflag  :  BOOLEAN; 

LoNerLiait,  OpperLialt  :  SINGLE); 


LABEL  localout; 

VAR  Integervalue  :  INTEGER; 

begin 


Val(NewVAlue, Integervalue, code) ; 

NewValuea  ; 

IF  code  <>  0  THOI 
begin 

MeWValueE  :>  errorl; 

Errorflag  FALSE; 

GOTO  localout; 
eitd; 

IF  (integervalue  <  lowerlialt)  OR  (integervalue  >  upperliait)  THEN 
begin 

Errorflag  :»  FALSE: 

NewValue2  :=  error2; 

GOTO  localout; 
end; 

localout  : 
end; 


<  This  procedure  checks  to  make  sure  that  the  filename  conforms  to  PC- DOS  standards 


procedure  CheckFlleNaae(VAR  naae,errorBessage:STRING); 

LABEL  localout; 

begin  {CheckFileNaae} 

errorflag  :«  TRUE; 
erroraessage  :■  ”; 


IF  length(naae)>12  Then 
begin 

Errorflag  :»  FALSE; 

Erroraessage  :*  'Length  >  12  Characters!*; 
GOTO  localout; 
end; 


IF  ((POSC  . '  ,naae)  >  9)  or  (POS(’ .  ’ ,naae)  * 
begin 

Errorflag  :«  FALSE; 

Erroraessage  :«  'Naae  >  8  Characters!'; 
GOTO  localout; 
end; 


0))  and  (lengthfnaae)  >*  9)  THEN 


IF  (Length(naNe)  -  P0S( ’ .  ’ ,naae))  >  3  THEN 
bMin 

Errorflag  :*  FALSE; 

Erroraessage  :«  'Extension  >  8  Characters!’; 
GOTO  localout; 
end; 


localout  :; 


end;  (CheckFileNaae) 


> 
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{ 


This  procedure  reatk  the  input  parameten  for  the  ntn. 


procedure  TMyApp.RunParaaeters; 
var 

Thoaas:  PView: 

Dialog:  POeaoOialog; 

R:  TRect; 


C:  Vtord; 

erroraessage , 

teapstrlng  :  STRING; 

yf irat, 

xleft, 

xright, 

xleft2, 

xrlght2, 

xleftS, 

xrlghtS  ;  BYTE; 
lower value  :  INTEGER; 


procedure  CheckRunParaaeters ; 

{  IhU  it  the  hxtal  ptooethire  that  cfaeckt  the  run  parameter  data  eniriet.  It  utet  ChceUnlcferOeiB  and 

ChtckRcalDala  with  the  appropriate  vahiet. 


VAR 

newcode, 

■InieuaM, 

■axieuNM  :  INTEGER: 
UpperBound  :  SINGLE; 

LABEL  Done; 


begin 

(  Fast  set  the  enotflag 

Runinput.erroreesaage  ; 
errorflag  TRUE; 

{  Then  check  each  of  (he  record  vahiei  suniag  with  (he  minimum  aumbtr  of  RDCciraJ  tennis  Minm.  If  an  error  occurs 

in  any  of  the  checks,  the  error  Qag  it  set  to  FALSE  before  exiting  the  procedure. 

CheckIntegerOat  a  ( Run  Input .  Nine ,  Runinput .  erroraessage , 

'Mlnlaue  Spectral  Tere  Not  Integer!’, 

’Mlniaua  Spectral  Tere  Out  Of  Range! * .errorflag,  3,  48); 

IF  errorflag  »  FALSE  THEN  GOTO  Done; 

VAL(RunInput. Mine, UlnleueH, code) ; 

CheckIntegerOataf  Runinput .  Maxe ,  Runinput .  erroraessage , 

’Haxlaua  Spectral  Tera  Not  Integer!', 

’Uaxlaua  Spectral  Tera  Out  Of  Range! ’ .errorflag,  HinlauaU,  51); 

IF  errorflag  ■  FALSE  THEN  GOTO  Done; 

CheckI  ntegerOata  (RunI  nput .  Npolnt  a ,  RunI  nput .  erroraessage , 

’Nuaber  Of  On-Axla  Points  Not  An  Integer!’, 

’Nuaber  Of  On-Axls  Points  Out  Of  Rangel ’ ,errorf lag,  99,  401); 

IF  errorflag  •  FALSE  THEN  GOTO  Done; 

CheckIntegerOata(RunInput  .NOpoints ,  Runinput .  erroraessage, 

'Nuaber  Of  Off-Axis  Points  Not  An  Integerl’, 

'Nuaber  Of  Off-Axis  Points  Out  Of  Range! ’ ,errorf lag,  0,  41); 

IF  errorflag  -  FALSE  THEN  GOTO  Done; 

CheckRealOatafRunInput  .DeltaR, Runinput  .erroraessage, 

'Radial  Step  Size  Not  A  Nuaber!’, 

'Radial  Step  Size  Out  Of  Range! errorflag,  0.001,  0.0201); 

IF  errorflag  >  FALSE  THEN  GOTO  Done; 

CbeckRealData(  RunI  nput  .Zaln, Runinput  .erroraessage, 

'Mlniaua  On-Axis  Output  Not  A  Nuaber!  ’ , 

'Mlniaua  On-Axis  Output  Of  Range! ' .errorflag,  0.001,  0.99); 

IF  errorflag  >  FALSE  THEN  GOTO  Done; 
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CheckRealData  ( Runinput .  Zaax ,  RunI  nput .  erroraessage , 

'Maxiaua  On-Axia  Output  Not  A  Nuaber!*, 

'Maxiaua  On -Axis  Output  Of  Range! ' ,errorf lag,  1.05,  2.00); 
IF  errorflag  =  FALSE  THEN  GOTO  Done; 


CheckRealData  ( Runinput .  Rgaax ,  Runinput .  erroraessage , 

'UaxlBua  Off-Axis  Range  Not  A  Nuaberl*, 

'Maxiaua  Off-Axis  Range  Out  Of  Range! ‘ ,errorf lag,  1.1,  10.0); 
IF  errorflag  »  FALSE  THEN  GOTO  Done; 


CheckI  ntegerOata  ( RunI  nput .  RadlusMax ,  Run  I  nput .  erroraessage , 
'Maxiaua  Off-Axis  Points  Not  A  Nuaber!', 

'Maxiaua  Off-Axis  Points  Out  Of  Rangel  ’  ,errorf lag,  10,  200); 
IF  errorflag  >  FALSE  THEN  GOTO  Done; 
done: 
end; 


begin 


gin 
Xleft 
Xrlght 
Xlefta 
Xrlght2 
Yflrst 


xieft  ♦  40; 
Xrlght  ♦  1 
Xleft2  ♦ 

4 


8; 


R.Asslgn(15,  75,  )rflrst  *  16); 
TEapStrlng  :>  'Maxiaua  i  Spectra 
Dialog  :<•  NewlPOeaoOlalog,  lnlt(R, 
with  Dialog*  oo 
be^ln 


(Min  #  to  50):’; 
' RunParaaeters ' ) ) ; 


(4  to  49):’,Thoaas))); 


n.Assign(xleft.  2,  xleft  ■*’  48,  3); 

Insert (New(PLaDel,  Inlt(R,  'Use  Tab  key  or  aouse  to  select  data  entry  field’, 
R.Assign(xleft2,  yflrst,  xrlght2,  yflrst  +1); 

Thoaas  :-  New(PInputLlne,  InIt(R,128)) ; 

Insert(Thoaas) ; 

R.As8ign(Xleft.  yflrst,  Xrlght,  yflrst  +1); 

Insert (NewIPlaoel, Inlt (R, ‘MinlauB  #  Spectra 
R.Assign(xleft2,  yflrst  ♦  i,  xrlght2,  yflrst  +  2); 

Thoaas  :=  New(PInputLine,  Init(R,l28)); 

Insert (Thoaas); 

R.Asslgn(xleft,  yflrst  ♦  1,  xrlght,  yflrst  ♦  2); 

Insert  (New(Plabel,Inlt(R,TeapStrIng,Trioaas))); 

R.Asslgn(xleft2,  yflrst  ♦  2,  xright2,  yflrst  ♦  3); 

Thoaas  :«  New(PInputLine,  Inlt(R,l26)); 

Insert (Thoaas); 

R.Assign(xleft,  yflrst  2,  xrlght,  yflrst  ♦  3); 

Insert(New(Plabel,Init(R, 'Nuaber  On-Axls  Points 
R.Assign(xleft2,  yflrst  +  3,  xright2,  yflrst  ♦  4); 

Thoaas  :=  New(PInputLlne,  Inlt(R,l28)); 

Insert (Thoaas); 

R. . 


(100  to  400):’, Thoaas))); 


R.Asslgn(xleft,  yflrst  ♦  3,  xrlght,  yflrst  +  4); 
Insert(New(Plabel,Inlt(R, 'Nuaber  Off-Axis  Points 
R. Assign (xleftE,  yflrst  ♦  4,  xright2,  yflrst  +  5); 
New(PinputLlne,  Inlt(R,l28)); 


(1  to  40):’, Thoaas))); 


5); 

(0.02  to  0.001):’, Thoaas))); 
♦  6); 


to  0.99):’, Thoaas))); 


R 

Thoaas 

Insert(Thoaasi; 

R. Assign (xleft,  yflrst  ♦  4,  xrlght,  yflrst 
Insert(New(Plabel,Inlt(R, 'Radial  Step  Size 
R.Asslgn(xieft2,  yflrst  >  5,  xright2,  yflrst 
Thoaas  :*  New(PlnputLine,  lnlt(R,l28)); 

Insert(Thoaas); 

R.A8slgn(xlefi.  )rflrst  *  5,  xrlght,  yflrst 
lnsert(New(Plaoel,Init(R, ’Min  On-Axis  Output 
R.Assign(xieft2,  yflrst  6,  xright2,  yflrst 
Thoaas  :«  New(PInputLine,  Init(R,128)); 

Insert (Thoaas); 

R.Assign(xlefi,  yflrst  +  6,  xrlght,  yflrst  +  7); 

Insert(New(Plabel,Inlt(R, 'Max  On-Axls  Output  (1.05  to  2.0):' ,ThoBas))); 
R.Assi(jn(xleft2,  yflrst  ♦  7,  xrlght2,  yflrst  8); 

Thoaas  :-  New(PlnputLlne,  Init(R,l28)); 

Insert (Thoaas); 

R.Asiign(xleft,  yflrst  ♦  7,  xrlght,  yflrst  +  8); 

Insert(New(PlaDel,Init(R, 'Max  Off-Axis  Range  (1.1 

R. Assign (xleft2,  yflrst  ♦  8,  xright2,  yflrst  ♦  9); 

Thoaas  New(PinputLine,  Inlt(R,l28)); 


f  6); 

(0.001 

♦  7); 


to  10.0):’, Thoaas))); 


Insert (Thoaas); 

R. Assign (xleft,  yflrst  ♦  8,  xrlght,  yflrst  ♦  9); 
Insert(New(PlaDel,Inlt(R, 'Max  Off-Axis  Points 
R.Assign(16,  yflrst  ♦  10,  56,  yflrst  +  11); 

Thoaas  :=  New(PlnputLine,  lnit(R,l28)); 


(10  to  200) :’, Thoaas))) ; 


Thoaas ) ) ) ; 
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Insert (Thoaas): 

R.Assian(l,  yflrst  ♦  10,  15,  yflrst  ♦  11); 

Insert(Ne«r(Plabel,Init(R,  ’Error  Message’  ,fhoaas))); 

{  These  are  the  escape  butlons  from  the  boa.  To  use  the  values  enieied  eiecutc  emOK.  To  discaid  the  changes 

execute  cmCanoel.  ) 

R.Aaslgn(xleft  ♦  g,  ytlrst  ♦  12,  xleft  ♦  19,  yfirst  ♦  14); 

Insert (NewtPButton,  Init(R,  '~0~k',  caOK,  bfOeTaultH): 

R.  Assign  (xleft  ♦  34,  yflrst  ♦  12  xleft  ♦  44,  yflrst  ♦  14); 

Insert (New(PButton,  Inlt(R,  'Cancel',  caCancel,  bfNorsal))); 

end; 

REP^T 

Olaloa* .SetOata(RunInput ) : 

C  :>  DeskTop*. ExecVlewf Dialog): 

If  C  o  caCancel  then  Dialog*. GetData(Runlnput): 

{  Check  to  see  whether  there  are  any  enoB.  If  an  error  has  oocuried,  then  do  not  exit  the  program.  ) 

CheckRunParaaeters ; 

UNTIL  ((C  *  caCancel)  or  (C  >  caOK))  and  (errorflag  »  TRUE); 

Dlspose(Dlalog,  Done); 
end; 


{ 


Hiis  pnxedure  is  used  to  input  or  change  transducer  related  parameters,  sudb  as  the  radius  of  curvature,  diameter, 
initial  oo-aiis  potwer,  and  frequency. 


) 


procedure  TttyApp.NewTransducer; 
var 

Thoaa82:  PView; 

01alog2:  POeBoOialog; 

R:  TRect; 

C;  Word; 

erroraesaage  :  STRING; 

xleft. 

xrlght , 

xleft2, 

xrlght2 , 

yflrat  :  BYTE; 


{  Ibis  is  the  local  procedure  that  checks  the  transducer  data  entries.  It  uses  CheeklnlcterDala  and  ChrehRealData 

with  the  appropriate  vakiea  ) 

Procedure  CheckTranaducer ; 

LABEL  Done; 

begin  {CheckTranaducer) 

{  Fust  set  the  error  message  and  error  flag.  Then  diecfc  each  of  the  data  catries  starting  with  the  diameter.  If  an 

error  occurs,  set  the  error  flag  to  FALSE  and  exit.  ) 

transducer.erroraessage  ”; 
errorflag  :=  TRUE; 

CheckReal0ata( transducer .dlaaeter .transducer .erroraessage , 

'Oiaaeter  Entry  Is  Not  A  Nuaber! ’ , 

’Olaaeter  Entry  Out  Of  Range! ' .errorflag,  0,  1); 

IF  errorflag  =  FALSE  THEN  GOTO  Done; 

CheckRealOata(transducer .RadiusOf C .transducer . erroraessage , 

’Focal  Length  la  Not  A  Nuaberl’, 

'Focal  Length  Out  Of  Range! ’ .errorflag,  0,  1); 

IF  errorflag  ■  FALSE  THEN  GOTO  Done; 

CheckflealOat a ( t  r ansducer . f  r equency , t  ransducer . erroraessage , 

'Frequency  Is  Not  A  Nuaber!', 

'Frequency  Out  of  Range! ' .errorflag,  0,  10); 

IF  errorflag  =  FALSE  THEN  GOTO  Done; 

CheckRealDatattransducer.pressure, transducer. erroraessage, 

'Pressure  Is  Not  A  Nuaberl’, 

'Pressure  Out  Of  Range! ' .errorflag,  0,  2000); 

done  : 

end;  {CheckTransducer} 
begin 

Yfirst  :*  2; 

Xleft  2; 

Xright  :«  Xleft  ♦  34; 

Xleft2  :*  Xright  *  1; 

Xright2  :=  Xleft2  ♦  10; 

R.A8Sign(5,  6,  58,  20); 

Dialogz  :=  Hew(POenoOlalog,  Init(R,  'Transducer  Characteristics’ )); 

with  Dialog2''  do 

begin 

R.Assign(Xleft.  Yfirst,  50,  Yfirst  ♦1); 

Insert (NewfPLaoel,  Init(R,  'Use  Tab  key  or  aouse  to  select  data  entry  field',  Thoaas2))); 
R.Assign(Xieft2,  Yfirst  +  2,  Xright2,  Yfirst  ♦  3); 

Thoaasz  :»  New(PInputLine,  Init(R,128)); 

Insert (Thoaas2); 

R.Assign(Xleft,  Yfirst  ♦  2,  Xright,  Yfirst  ♦  3); 

In8ert(New(Plabel,Init(R, 'Diaaeter  (Meters:  0  ..  i ) : ’ ,Thoaas2))); 

R.Assign(xieft2,  Yfirst  +  3,  Xright2,  Yfirst  ♦  4); 

Thoaas2  :«  New(PlnputLine,  lnit(R,128)); 

Insert (Thoaas2) ; 

R. Assign (xleft,  yfirst  ♦  3,  xright,  Yfirst  +  4); 

Insert(New(Plabel,Inlt(R, 'Focal  Length  (Meters:  0  ..  1 ) : ’ ,Thonas2))); 

R.A6sign(xleft2,  Yfirst  +  4,  Xrightz,  Yfirst  ♦  5); 
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ThoBas2  :*  New(PInputLlne,  Inlt(R,128)); 

Insart (Thoaas2 ) : 

R.Assian(xleft,  yflrst  ■M,  xright,  Yflrat  ♦  5); 

Inaert(Ne«i(Plabel,Inlt(R, 'initial  Frequency  (IMz:  0  ..  10):',Thoaas2))); 
R.Aasign(xleft2,  Yflrat  ♦  5.  Xright2,  Yfirat  ♦  6); 

ThoaasS  :*  New(PInputLlne,  Init(R,128)); 

Insert  (TbMaa2); 

R.Aaslgn(xleft.  yflrat  ♦  S,  xright,  Yflrat  *  6); 

Inaert(New(Plwel,Init(R, 'Initial  Pressure  (kPa:  0  ..  2000):’,Th<Mias2))): 
R.Assign(16, Yflrat  ♦  7,  49,  Yflrat  ♦  8); 

ThoaasS  :>  New(PinputLine,  Init(R,128)); 

Insert  (Tho)ias2); 

R.Aasign(1 , Yflrat  ♦  7,  15, Yflrat  ♦  8); 

Insert(New(Plabel,Inlt(R, 'Error  Message’ ,Thoeas2))): 


(  Theie  are  the  eacape  buttons  from  the  bca.  To  use  the  vakiet  entered  execute  canOiL  To  discard  the  changes 

esecutc  cesCaaoel. 


> 


R.Aasign(10,  Yflrat  ♦  9,  20,  Yflrat  11); 
Inaert(New(rautton,  lnlt(R,  "O'k*,  caOK,  bfOef ault ) ) ) ; 
R.Aasign(35,  Yflrat  9,  45,  Yflrat  ♦  11); 
lnsert(New(reutton,  Inlt(R,  'Cancel*,  ceCmcel,  bf Norsal ) ) ) ; 


end; 

REPEAT 

01alog2* .SetOata( Transducer) ; 

C  :«  DeskTop*. ExecView(0ialog2); 

If  C  <>  caCancel  then  Dlalog2*.GetOata(Transducer); 

(  Check  to  tee  whether  there  are  any  ertos.  If  an  error  bat  occurred,  then  do  not  exit  the  program.  ) 

CheckTransducer; 

UNTIL  (((C  s  caCancel)  or  (C  =  caOK))  and  (errorflag  »  TRUE)); 

0ispose(Dialog2 ,Done) ; 
end; 
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(  This  procedure  is  used  lo  select  the  type  and  souice  of  the  ioitia]  pressure  disuibutioo  at  the  face  of  the  transducer. 

The  options  allow  the  user  to  enter  the  data  or  select  a  file. 

procedure  TMyApp.SelectBoundary ; 
var 

ThoBas2:  PVlew: 

0ialog2:  POeaoOialoa; 

R:  TRect; 


C:  Word; 
erroraesaage 
xleft, 
xrlght , 
xleft2, 
xrlghti , 

Xlef tBi , 
YleftBI , 
XRlghtBI , 
YRightBl , 
XleftB2. 
XRightB2, 
YLeTtB2, 
YRightB2, 
yflrat 
textflle 


STRING; 


BYTE; 

TEXT; 


begin 

Yfiret  :*  2; 

Xleft  :*  2; 

XleftBI  6; 

YleftBI  :»  5; 

XRightBI  XleftBI  >14; 

YRightBl  :=  YleftBI  ♦  3; 

XleftB2  XRightBI  *  6; 

YleftB2  :=  YleftBI; 

XRightB2  XLeftBi  >16; 

YRightB2  :»  YleftB2  >  2; 

Xright  ;=  Xleft  >  34; 

Xleft2  ;»  Xright  >  1; 

Xright2  :*  Xleft2  >10; 

R.Assign(S,  6,  55,  161; 

Oialogz  :*  NeeffPOeaoOialog,  Init(R,  'Transducer  Pressure  Distribution  Options’)); 

with  0ialog2*  do 

begin 

R.Assign(Xleft,  Yfirst,  46,  Yfirst  >  i); 

Insert (NewfPLabel,  Inlt(R,  'To  Select  Option,  Press  Highlighted  Letter’,  Thoaas2))); 

{  This  set  of  radio  bottoas  allow  the  selection  of  the  values  used  in  a  previous  run,  a  uniform  disirfcution,  or  a  set  of 

dta  values  contained  in  a  file. 


R.AssignfXleftBI .YleftBI .XRightBI , YRightBl ) ; 

Thoeasz  :«  New(PRadioButtons,Inlt(R, 

NewSIteaf ’*P"revious’ , 

NewSIteaf  ’  K’nlfora  ’ , 

Ne«iSItea(’~N~ew’ , 
nil))) 

) )  f 

Insert (Thoaas2); 

R.Assign(XleftB1-4,  YleftBI  -  1,  XRight;r>‘>4,  YLeftBI); 
Insert(New(PLabel,  Inlt(R,  ’Pressure  Cistri^ution’ ,  Thoaas2))); 

This  allows  the  selection  of  a  previous  file  or  the  ^  .-'ation  of  new  file. 

R .Assign (Xlef tB2,YleftB2. XRlghtB2 , YRlghtB2) ; 

Thooasz  :>  New(PRadio6uttons,Init(R, 

NewSItea(’Use  ~F“ile’, 

MewSIteaf ’"E“nter  Data’, 
nil)) 

)); 

Insert (Thonas2); 

R.Assign(XleftB2  >  4,  YleftB2-  1,  XRightB2,  YleftB2); 
Insert(New(PLabel,  Init(R,  ’Source’,  ThoaasZ))); 
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{  These  are  the  escape  buttons  from  the  hex.  To  use  the  values  entered  esecute  emOK.  To  discard  the  changes 

eiecute  cmCanoel.  ) 

R.Asslgn(7,  Yflrst  ♦  7,  17,  Yfirst  ♦  9); 

Insert(New(PButtoa,  Init(R,  caOK,  bfDefault))); 

R.A88lgn(29,  Yfirat  ♦  7,  39,  Yfir*t  ♦  9); 

In8ert(New(rautton,  Inlt(R,  ‘Cancel’,  caCancel,  bfNoraal))); 

end: 

REPEAT 

Dlalog2* .SetData(8oundaryChoice) ; 

C  DeskTop*. ExecVlew(Dlalog2): 

If  C  <>  caCancel  THEN  Dialogz*.6etOata(BoundaryCholce): 

UNTIL  (((C  *  caCancel)  or  (C  «  caOK))); 

Dispose (Dialog2 , Done) ; 
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{  This  procedure  is  used  to  enter  and  check  the  name  assoicated  with  an  existii^  data  6le. 

procedure  TMyApp.EnterFlleNaae; 
var 

Thoaas2:  PView: 

0ialog2:  POeaoOialog: 

R:  TRect; 

C:  Word; 

erroraessage  :  STRING; 

xleft, 

xright , 

xleft2, 

xrightz, 

ButtonLength , 

OKXlocatlon, 

OKYlocation, 

CXlocatlon, 

CYlocatlon, 
y first  ;  BYTE; 
textfile  ;  TEXT: 

Dir Info  :  SearchRec; 

begin  {EnterFlleNaae} 

Xleft  :■  2; 

Xright  :«  Xleft  ♦  22; 

Yfirat  :»  2; 

Xleft2  :»  Xright  ♦  2; 

Xright2  Xleft2  *  30; 

ButtonLength  :=  10; 

OKXlocatlon  :*  10; 

CXlocatlon  35; 

OKYlocation  :=  6; 

CYlocatlon  OKYlocation; 


ButtonLength  :=  10; 

OKXlocatlon  :*  10; 

CXlocatlon  35; 

OKYlocation  :=  6; 

CYlocatlon  OKYlocation; 

R.Asslgn(5,  6,  65,  15); 

0ialog2  :=  New(P0eao0ialog,  Inlt(R,  'Pressure  Distribution  File  Naae’)); 

with  0131002"  do 

begin 

R.Assign(Xleft2,  Yflrst,  Xright2,  Yfirst+1); 

ThoaasS  :»  Ne«(PlnputLlne,  Inlt(R,128)); 

Insert (Thoaas2) ; 

R.Asslgn(Xleft,  Yflrst,  Xright,  Yflrst  +1): 

Insert(New(Plabel,Init(fl,  "F'ile  Naae  (<40  Char) ’ ,Thoaas2) ) ) ; 

R.Assign(Xleft2,  Yfirst+2,  Xright2,  Yfirst+S); 

Thoaas2  :«  New(PInputLlne,  Inlt(R,128)); 

Insert (Thoaa82 ) ; 

R.Asslgn(Xleft,  Yfir8t+2,  Xright,  Yflrst  3); 

Insert(New(Plabel,Init(R, 'Error  Message’ ,Thoaas2))); 

These  are  (be  escape  buttons  from  the  bcK.  To  use  the  vahies  entered  ciecuie  cmOK.  To  discard  the  changes 
execute  cmCatKcl. 

R.  Assign  (OKXlocatlon,  OKYLocation,  OKXlocatlon  ButtonLength,  OKYlocation  *  2); 

Insert (New(PButton,  lnlt(R,  '~0~k',  caOK,  bfOefault))); 

R.Asslgn(CXlocation,  CYlocatlon,  CXlocatlon  ButtonLength,  CYlocatlon  *  2); 
lnsert(New(PButton,  lnlt(R,  'Cancel',  caCancel,  bfNoraal))); 

end; 

REPEAT 

01alog2"  .SetOata(FileNaBe) : 

C  :«  DeskTop*. ExecView(0ialog2 ) ; 

if  C  <>  caCancel  then 

begin 

skipflag  :=  FALSE: 

0ialog2' . GetOata ( FlleNaae ) ; 
checkf  llenaae  ( f  ilenaae  .naae,filenaBe.  erroraessage ) : 

IF  (C  >  caOK)  and  (FileNaae.Naae  «  ”)  THEN  errorflag  :=  FALSE 

ELSE 

IF  Errorflag  <>  FALSE  THEN 
begin 

FindFirst  (f  ilenaae .  naae ,  AnyFlle ,  Oir  Inf  o ) ; 

IF  DosError  <>  0  THEN 
begin 

errorflag  :»  FALSE; 

f ilenaae. erroraessage  :=  'File  ’-•'f ilenaae. naae-t-  '  Not  Found!’; 
end; 
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else 


end; 
end 
begin 

errorflag  :=  TRUE; 
sklpflag  TRUE; 
end; 

UNTIL  (((C  s  ceCancel)  or  (C  >  caOK))  and  (errorflag 
IF  C  ■  caCancel  THEN  BoundaryCtiolce.RadioButtonOataA 
01spose(01alog2,Done) ; 

end;  {EnterFlleNaie} 


a  TRUE)) 
:»  0; 
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{ 


This  procedure  is  ittcd  to  enter  the  name  of  a  data  fik  sriucb  is  to  be  created,  the  number  of  data  points  which  will 
be  contained  in  the  file,  and  check  the  file  name  to  make  sure  it  oonfonns  oi  DOS  standarik. 


) 


procedure  TllyApp.CnterFileNaae2; 
var 

Thoaas2:  PView; 

01alog2:  POeaoOialog: 

H:  TRect; 

C:  Word; 

erroraessage  :  STRING; 

xleft, 

xright, 

xleTt2, 

xrlghtz, 

ButtonLength , 

OKXlocatlon, 

OKYlocatlon, 

CXlocatlon, 

CYlocatlon , 
yf Irat  ;  BYTE ; 

textflle  :  TEXT; 

Procedure  CheckNuaberPoints; 

(  This  piocedure  checks  the  number  of  points  to  make  sure  it  is  an  integer  and  witlun  range.  } 


LABEL  Done; 

begin  {CheckNuaberPoints} 

fllenaBe2. erroraessage 
errorflag  TRUE: 

CheckInteger0ata(fllenaae2.NuaberPolnts,filenaae2.erroraessage, 
'Nuaber  Of  Points  Must  Be  >  O’, 

’Nuaber  Of  Points  Out  Of  Range! ’ .errorflag,  0,  100); 
end;  {CheckNuaberPoints} 


begin  {EnterFlleNanea} 

Xleft  :»  2; 

Xright  :»  Xleft  22; 

Yfirst  :»  2; 

Xleft2  :»  Xright  ♦  2; 

Xrlght2  :»  Xleft2  *  30; 

ButtonLength  10; 

OKXlocatlon  i0; 

CXlocatlon  35; 

OKYlocatlon  7; 

CYlocatlon  :=  OKYlocatlon; 
R.Asslgn(5,  6,  65,  16); 

0ialog2  New(P0eao0ialog,  Inlt(R, 
with  Dlalog2*  do 


’New  Pressure  Distribution  File  Naae’)); 


begin 

R. Assign (Xleft2,  Yfirst,  Xrlght2,  Yflrst+1); 

ThoaasB  :*  New(PlnputLine,  Init(R,f28)); 

Insert (Thoaas2); 

R. Assign (Xleft,  Yfirst,  Xright,  Yfirst  *1); 
Insert(New(Plabel,Init(R,’-F~iie  Naae  (<40  Char)’,Thoaas2))); 
R.Assign(Xleft2,  Yflrst-M,  Xright2,  Yflr8te2); 

Thoaas2  :*  New(PInputLine,  Init(R,128)); 

Insert(ThoBas2); 

R.AsslgnfXleft,  Vflrst+1,  Xright,  Yfirst  ♦  2); 
Insert(New(Plabel,Init(R, ’~N~uaber  of  points’, Thoaas2))); 


R.Assign(Xleft2,  Yfirst+3,  Xright2,  Yfirst+4); 
Thoaasz  :>  New(PInputLine,  Init(R,128)); 
Insert(Thoaas2); 

R.Assign(Xleft,  Yfirst+3,  Xright,  Yfirst  ♦  4); 
Insert(New(Plabel,Init{R, 'Error  Message' ,Thoaas2))); 


{  These  are  the  escape  buttons  from  the  bca.  To  use  the  values  entered  execute  emOK.  To  discard  the  changes 

execute  cmCanoel. 


) 


R.A88ign(0KXlocation,  OKYLocation,  OKXlocation  *  ButtonLength,  OKYlocatlon  -i-  2); 
Insert(New(PButton,  Init(R,  ’~0~k’,  caOK,  bf Default ))) ; 

R.Asslgn(CXlocation,  CYlocatlon,  CXlocatlon  ButtonLength,  CYlocatlon  *  2); 
Insert (New(PButton,  Inlt(R,  ’Cancel’,  caCancel,  bfNoraal))); 
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•nd; 

REPEAT 

skipflag  :3  FALSE: 

01alog2* .SetOata(FlleNaBe2) ; 

C  :=  DeskTop* .ExecView(Dlalog2): 

If  C  <>  caCancel  then 
begin 

Dialoa2* .GetDatafFileNaaeS ) ; 

checkfl lenaae (f ilenaae2 . naae , f ilenaae2 . erroraessage ) ; 


IF  ((C  «  caOK)  and  (FileNaae2.Naae  >  ”))  THEN 
begin 

errorflag  FALSE; 

filenaaez.erroraessage  :>  'No  File  Naae  Entered!’; 
end 

ELSE 

IF  errorflag  ■  TRUE  THEN  CheckNuaberPoints; 
end  else 

begin 

errorflag  :>  TRUE; 
skipflag  :«  true; 
end; 

UNTIL  ((C  *  caCancel)  or  ((C  -  caOK)  and  (errorflag  ^  TRUE))) 
IF  C  >  caCancel  THBI  BounoaryChoice.RadioButtonOataA  0; 
Dispose(0ialog2 ,0one) ; 

end;  {EnterFlleNaae2) 
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{  Hiu  procedure  ii  OKd  for  the  off-axis  data  entry. 


procedure  TUyApp.OffAxisOata; 
var 

Thoaas2:  PVlew: 

01alog2:  POeaoOlalog; 

R:  TRect; 

C:  Mord; 
teapwal, 

erroraessage  :  STRING; 
xlef  t , 

xrlght2, 

yflrat  :  BYTE; 

{  This  it  the  local  procedure  that  checks  the  transducer  data  entries.  It 

with  the  appropriate  vahies. 


ChecklniefcrDnia  and  ChrcUtealDaU 


Procedure  CheckOffAxlsOata; 

VAR 

Binval, 

aaxval  ‘.SINGLE; 

LA8Q.  Done; 

begin  {CheckOffAxlsOata} 

(  Fvst  set  the  error  message  and  error  Then  check  each  dau  entry.  Snce  this  procedure  is  caOed  each  time  the 

off  axis  index  changes.  If  an  error  occurv  set  the  error  flag  to  FALSE  and  exit. 


} 


) 


} 


ZU(offaxislndex|  .erroraessage  “; 
errorflag  :»  TRUE; 

VAL(ZU(offaxislndex  -1  ] . location, ainval, code) ; 

VALfRunlnput .Zaax, aaxval , code) ; 

CheckRealOata(ZU(offaxisindexi  .location, ZU(offaxisindex]  .erroraessage, 

'Off-Axis  Entry  Is  Mot  A  Nuaber!’,  'Previous  Entry 

ZUloffaxlsindex  -  1). location  *  ’  Entry  Out  Of  Range! errorflag,  ainval,  aaxval); 
IF  errorflag  *  FALSE  THEN  GOTO  Done; 
done  : 


end;  {CheckOffAxlsOata} 


begin 

Yfirst  ;*  2; 


lalog,  Init(R,  'Off  Axis  Location  teapval)); 


Xleft  :>  5; 

Xright  Xleft  *  38; 

Xlert2  :»  Xright  ♦  1; 

Xright2  :>  Xleft2  *  10; 

STRfoffaxislndex, teapval ) ; 

R.Assign(5,  6,  65,17); 

OialogS  :*  NewfPOeaoOia 
with  01alog2*  oo 
begin 

iTMsign (Xleft.  Yfirst,  Xleft  ♦  50,  Yfirst  ♦  1); 

Insert(New(PLabel,  Init(R,  'Use  Tab  key  or  aouse  to  select  data  entry  field’,  Thoaas2))); 
R.Assign(xleft2,  Yfirst  ♦  2,  XrightE,  Yfirst  ♦  3); 

Thoaasz  :«  NewfPInputLlne,  Inlt(R,128)); 

Insert (Thoaas2 ) ; 

R. Assign (Xleft.  Yfirst  ♦  2,  Xright,  Yfirst  ♦  3); 

lnsert(New(Plaoel,lnit(R, 'Output  Point  (>  previous  value  <  Runinput.Zaax  '«-’):’,Thoaas2))); 
R.Asslgn(16, Yfirst  ♦  4,50,  Yfirst  ♦  5); 

ThoeasS  :*  New(PinputLine,  Init(R,l28)); 

Insert (Thoaas2 ) ; 

R.Assign(1, Yfirst  ♦  4, 15, Yfirst  ♦  5); 

Insert(New(Plabel,Init(R, 'Error  llessage’ ,ThoBas2))); 

These  aie  the  escape  buttons  Croat  the  boa.  To  use  the  vahies  entered  execute  cmOK.  To  discard  the  changes 
execute  onCancel.  } 


R.AsslgnflO,  Yfirst  ♦  6,  20,  Yfirst  *  8); 

Insert (New(PButton,  Init(R,  ’~0~k’,  caOK,  bfDefault))); 
R.Assign(35,  Yfirst  ♦  6,  45,  Yfirst  ♦  8); 

Insert(New(rautton,  Init(R,  'Cancel',  caCancel,  bf Noraal ) ) ) ; 
end; 
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REPEAT 

01alog2'‘  .SetData(ZU(Off  Axlslndex) ) ; 

C  OeakTop* .ExecView(01alog2) : 

if  C  <>  caCancel  then  01alog2*.Get0ata(ZU(0ffAxisIndex) ) ; 

{  Check  to  see  whether  there  are  any  enon.  If  ao  error  has  occurred,  then  do  oot  exit  the  program.  ) 

CheckOffAxlsOata; 

UNTIL  (((C  '  caCancel)  or  (C  *  caOK))  and  (errorflag  =  TRUE)); 

Dispose ( Dialog2 , Done ) ; 
end; 
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{  This  piooedure  chechs  each  location  and  pressure  pair  to  make  sure  they  are  sirithin  range. 


) 


procedure  THyApp.BoundaryAxisOata; 

VAR 

Thoflas2:  PView: 

01alog2:  POeaoOialog; 

R:  TRect; 

C:  WORD; 
teapval, 

erroraessage  :  STRING; 

xleft, 

xrlght , 

xlefta, 

xrlght2, 

yfirat  :  BYTE; 

procedure  CheckBoundaryAxlsOata; 

(  This  procedure  checks  the  normalued  pressure  or  the  radial  input  point  to  assure  that  the  nomialized  pressure 

less  than  IjO  and  that  the  radial  input  point  is  less  than  the  transducer's  transverse  radius. 


) 


VAR 

ainval, 

aaxval  tSINGLE; 

LABEL  Done; 

begin  {CheckBoundaryAxlsOata} 


BoundaryCholceValues{offaxlslndex|  .erroraessage  “; 
errorflag  :=  TRUE; 
ainval  :=  0; 


VAL (Transducer .Oiaaeter , aaxval , code ) ; 
aaxval  ;=  aaxval/2.0; 

CheckRealOata(BoundaryCholceValues(offaxlslndex|  .location, 
BoundaryCholceValues(offaxlslndex)  .erroraessage, 

•Oft-Axls  Location  Entry  Is  Hot  A  Nuaber! ' , 

'Location  '  >  BoundaryCnolceValuesloffaxlsindex]  .location-*’ 
’  Out  Of  Range!',  errorflag,  ainval,  aaxval); 

IF  errorflag  =  FALSE  THEN  GOTO  Done; 


ainval  :*  0; 
aaxval  l ; 

CheckRealOata(BoundaryCholceValues( of f  axlslndex]  .NoraalizedP, 
BoundaryCholceValues [of f axlslndex | . erroraessage , 

'Pressure  Is  Not  A  Nuaber! ' , 

'Pressure  Entry  Out  Of  Range!',  errorflag,  ainval,  aaxval); 

IF  errorflag  >  FALSE  THEN  GOTO  Done; 


done  : 

end;  {CheckBoundaryAxlsOata} 


begin 

Yflrst  :«  2; 

Xleft  8; 

Xrlght  Xleft  28; 

Xleft2  :■  Xrlght  *  1; 

Xright2  :*  Xleft2  *  10; 

STR(off axlslndex  -*  1, teapval); 

R.Asslgn(S,  6,  65,18): 

OlalogS  New(POeao6ialog,  Inlt(R,  'Pressure  Distribution  Location  •  '*  teapval)); 

WITH  01alog2*  00 
begin 

R.Asslgn(Xleft-3,  Yflrst,  Xleft  ♦  50,  Yflrst  *1); 

Insert (New(PLabei,  Inlt(R,  'To  Change  The  Value,  Press  Highlighted  Letter',  Thoaas2))); 
R.Asslgn(Xleft2,  Yflrst  2,  Xrlght2,  Yflrst  *  3); 

Thoaas2  :>  Ne«(PInputLlne,  inlt(R,l28)); 
lnsert(Thoaas2); 

R.Asslgn(Xleft,  Yflrst  +  2,  Xrlght,  Yflrst  ♦  3); 

Insert(New(PlaDei,Inlt(R, 'Radial  ~P~olnt  (B)’,thonas2))); 

R.Asslgn(Xleft2,  Yflrst  ♦  3,  Xrlght2,  Yflrst  ♦  4); 

ThoaasS  New(PlnputLlne,  lnlt(R,l28)); 

Insert (Thoaas2); 
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R.Asslgn(Xleft,  Yfirst  ♦  3,  Xright,  Yfirst  ♦  4); 

Insert(New(Plabel,Init(R, ’’N'oraalized  Pressure’ ,Thoaas2) ) ) ; 

R.Assign(l6,Yflrst  ♦  5,59,  Yfirst  ♦  6); 

Th(Mias2  New(PlnputLine,  Inlt(R,  128) ) ; 
lPsert(Thoaas2) ; 

R.Asslgn(l , Yfirst  ♦  5, 15, Yfirst  ♦  6); 

Insert(Nei«(Plabel,Init(R,  'Error  Message’  ,Thoaas2))): 

(  These  aic  the  escape  buttons  Grooi  the  bat.  To  use  the  values  entered  caecuie  cmOK.  To  discard  the  chaopcs 

eseaite  cmCancel.  } 

R.Asslgn(10.  Yfirst  ♦  7.  20,  Yfirst  ♦  9); 
lnsert(Nei«(rautton,  Init(R,  ’~0~k’,  caOK,  bfDefauIt))); 

R.Assign(35,  Yfirst  ♦  7,  45,  Yfirst  ♦  9); 

Insert(New(PButton,  Init(R,  ’Cancel*,  caCancel,  bfNoraal))); 

end: 

REPEAT 

0181002* .SetOata(BoundaryCholceValues(Of fAxlslndex] ) ; 

C  :*  DeskTop*. ExecVlew(Dialog2); 

IF  C  <>  ceCancel  THEN 

0i8log2*  .GetOata(BoundaryCholceValues(Of f Axisindex] ) ; 

CheckBoundar yAxisOata : 

UNTIL  (^C  «  caCancel)  OR  (C  *  caOK))  ANO  (errorflag  ■  TRUE)); 

0ispose(0ialog2,0one); 


a 
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procedure  TMyApp.NewWindow; 
ver 

Window:  POeaoWindow; 

R:  TRect; 
beoln 

inc(WlnCount); 

R.Assign(0,  6,  45,  13); 

R.liove(Randow(34),  Randoa(ll)); 

Window  ;»  New(POewoWindow,  Init(R,  'Oeao  Window*,  WinCount)); 
DeskTop* .  Insert  (Window) ; 
end; 

var 

MyApp:  TllyApp; 

{  This  procedure  b  used  to  initialize  the  radio  button  choices  and  the  file  names. 


procedure  InitiateBoundary_FileNaae; 


in 


ITH  BoundaryCttolce  DO 
begin 

RadioSuttonOataA  :>  0; 
RadloButtonOataB  :■  0; 
end; 

WITH  FlleNaae  00 
begin 

Naae  ; 
ErrorMessage  ;= 
end; 

WITH  FileNane2  DO 
begin 
Naae 

NuaberPoints 
ErrorUessage 
end; 
end; 
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{  This  procedure  saves  the  parameter  values  selected  for  use  as  defaults  lor  the  nest  luo. 

procedure  WriteOefault; 

VAR 

textfile  :  TEXT; 


begin  fHrlteOefault} 
Asalgnftextf lie, ’Default .dat ' ) ; 
Rewrlte(textf lie) ; 

With  Transducer  do 


begin 

Writelnftextf ile,Dlaaeter) ; 
Writelnftextf lle,RadlusOfC) ; 
Writelnftextf lie, frequency) ; 
Writeln(textf lie, pressure) ; 
end: 

With  Runinput  Do 
begin 

Writelnftextf  lle,HlMt) ; 
Writelnftextf  lle,ltAXM) ; 
Writeln(textfile,Npolnte) ; 
Writelnftextf lie, NOpoints); 
Writelnftextf lie, OeltaR); 
Writelnftextfile,2ein); 
Writelnftextflle,Zaax) ; 
Writelnftextflle,Rgaax) ; 
Writelnftextf  lie, RadlusUax) ; 
end; 

For  n  0  to  50  do 
WITH  ZUfn]  Do 
begin 

Writelnftextfile,locatlon) ; 


end; 

Writelnftextf ile,NuaberBoundaryPoints ) ; 

FOR  n  :s  0  TO  NuaberBoundaryPolnts  -  1  00 
begin 

VALf8oundaryCholceValues(nI  .location,teap1  ,code); 

VAL  f  BoundaryChoiceValues  { n  j .  NoraalizedP ,  t  eip2 ,  code ) ; 
Writelnftextflle,Teep1 :12;8) ; 

Writelnftextf ile,Te«p2; 12:8) ; 
end; 

Closeftextfile) ; 
end;  {WriteOefault} 
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{  This  procedure  rearfc  the  parameter  values  selected  for  use  as  defaults  for  the  next  run. 

procedure  ReadOefault; 

VAR 

textfile  :  TEXT; 


begin  (ReadOefault) 
Asslgnftextflle, ’Oefault.dat’ ) ; 
Reset(textflle) ; 

With  Transducer  do 
begin 

Readln(textflle,Dlaaeter) ; 
Readln(textflle,RadlusOfC) ; 
Readln(textf lie, frequency) ; 
Readln(textf lie, pressure) ; 
erroraessage 
end: 

With  Runinput  Oo 
begin 

Readlnftextf ile.MINM) ; 
Readlnftextflle.MAXM) : 
Readlnftextf lie, Npointa); 
Readlnftextf lie, NOpolnts) ; 
Readlnftextf lie, DeltaR) ; 
Readlnftextf lie, Zaln) ; 
Readln(textflle,ZBax) ; 
Readln(textflle,Rgaax) ; 
Readlnftextf lie, RadlusMax) ; 
erroraessage  ‘ ; 
end; 

For  n  0  to  50  do 


WITH  ZU[n|  Oo 
begin 

Readlnftextfile, location) : 
erroraessage  := 
end: 

READLN ( textfile , NuaberBoundaryPolnt s ) ; 

FOR  n  0  TO  NuaberBoundaryPolnts  •  1  00 
begin 


tgln 

REAOLN  f  text f lie , BoundaryCholceValues ( n ] . location ) ; 
REAOLN  f textf i le , BoundaryCholceValues ( n ] ,  NoraallzedP ) ; 
BoundaryCholceValuesfnl .erroraessage  :»  **; 


end; 

FOR  n  :»  NuaberBoundaryPolnts  TO  100  00 
WITH  BoundaryCholceValuesfnl  00 
begin 

location  :*  ” ; 

NoraallzedP  := 


erroraessage 

end; 

Close(textfile) ; 
end;  (ReadOefault) 
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{  ThU  procedure  wrto  the  input  dau  file  (or  the  LFOCUS  run. 

procedure  WrltelnputData; 

VAR 

textflle  :  TEXT; 

teapstrlng  :  STRING; 

teflpval , 
radlusaaxval , 
deltarOut , 

radlusval  :  DOUBLE; 

teaplnt  :  INTEGER; 


begin  (WrltelnputOata) 

Aaalgn(  textflle,  ‘  In^.dat* ) ; 

Rewrlte(textf lie) ; 
vrltelnitextf lie, M  ’): 

Witeln  (textflle, '  1  .M-B’ ) : 

NTlteln  ( textflle ,  Runlnput .  RadlusUax ) ; 

Writeln  (textflle  ,RunInput .  f^iax^  ’  00  ‘ ) ; 

Nriteln(textf  lie,  *2.Sd-l4' ) ; 

Wrlteln(textf  lle,Tran8ducer  .pressure-^'O-^S’ ) ; 
«rlteln(textf  lie,  ’9.96d-»’2' ) ; 
llflteln(textf  lie,  ‘5.SOO* ) ; 

Writeln  (textflle ,  Tranaducer .  RadlusOf  ’  00 ' ) ; 

VAL  (Transducer .  Oiaaeter ,  teepval ,  code ) ; 
teepval  :>  teapval/2.0; 

STR  ( tenpval :  9 : 6 ,  teapstrlng ) ; 
Wrlteln(textflle,teapstrlng'V>oo* ) ; 

Wrlteln(textf  lie, Runlnput  .OeltaR-r^'oO* ) ; 

Wr iteln( textflle, Runlnput. Zaln-f* 00' ) ; 
Wrlteln(textflle,Runlnput.ZBax‘*^’00’ ) ; 
Writeln(textflle,Runlnput  .Npoints) ; 
Wrlteln(textfile,Runlnput.NOpoints) ; 

VAL  ( Runlnput .  NOpolnts  iTeaplnt ,  code ) ; 

For  n  1  to  teapint  Oo 

WITH  ZU(n]  Oo 

begin 

writeln(textflle,locatlon4’00* ) ; 
end; 

With  Runlnput  Oo 
begin 

Writeln(textf ile.HINH) ; 

Writeln(textf  lle.lMXH) ; 
end; 

Writeln(textf  lie,  ’  1 ,4920+3' ) ; 

Writeln  { textf  lie,  Transducer  .frequency-*’  ’  D-*€  ’ ) ; 

Writeln  ( text  f  lie  ,NuaberBoundaryPolnt8 ) ; 

FOR  n  0  TO  NuaberBoundary Points  •  1  00 
begin 

VAL(BoundaryCholceValues(n]  .locatlon,teap1  ,code); 
VAL  ( Boundar yCholceValues  { n  j .  Noraal  izedP ,  teap2 ,  code ) 
Writeln(textfile,Teap1  :t2;8, '  ’ ,Teap2;t2:8); 

end; 

Close (textflle); 
end;  {WMteInputOata} 
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{  This  procedure  wrhcs  •  file  which  ooauins  the  kxatioa  aonnalized  pieaore  pair  that  is  used  to  describe  the  initial 

amplitttde  tfistriNiliaa  at  the  face  of  the  transducer.  } 

procedure  WriteBouadaiyFile; 

VAR 

lenfile  :  TEXT; 
n  :  INTEGER; 

Tempi, 

Temp2  :  DOUBLE; 

begiii  (WfiteBoundaryFile} 

As«gn(teitfile,HleName2.name); 

Re«nite(tciifile); 

Wrtleln(lestfilejMleName2,NumberPDiiiu); 

FOR  n  :•  0  ID  NumberBoundaiyPoinis  -  1  DO 
begin 

VAL(BaundaryC]ioioeVahici(n].lacaltoii,tefflpl.oode); 

VAL(BooiidafyChoiceVahies(o|J<ormalaedP4emp2^Mde); 

Writelii(lenfile,Templ:12:8); 

Writein(teitfile,Tefflp2: 1 2:8); 
end; 

CLOSEfteitfik); 

end;  {WriteBcundaryFile} 

begin  (Menu) 

ZU(0|.localioo  ;=  T).0'; 

ReadDefailt; 

InitiateBouodary_FileNaine; 
errorflag  :»  TRUE; 

MyAppJnit; 

MyApp.Run; 

MyAppDooe; 

WriteDe  fault 

IF  BoundaryChoice.RadioBuitoaOataA  =  2  THEN 
If  BoundatyChoice.RadioBuctoaDataB  =  I  THEN  WriteBoundaiyFile: 

Write  InputElata; 
end.  (Menu) 
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APPENDIX  C 
LFOCUS.F 
PROGRAM  LISTING 
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Some  addiiioaal  terms  used  in  the  program 


ZMIN 

ZMAX 

RGMAX 

Rpnin 

Rmax 

IflRESHOLD 

MINM 

MAXM 

NPOINTS 

NDFOINTS 

Number  BP 

TPOOm 

lOUTMAX 

MARRAYINDEX 

IMAX 

IMAXl 

IMAX4 

MAXMIN 

MAXMAX 

RESTART 

K 

L 

COUNTER 

CASE 

NSTEP 

Pi 

WOPi 

PO 

Ctaega 
BovetA 
BcverM 
Number  RP 
NcwDeluR 


-  Minimum  Z  value  for  output  purposes 

-  Maximum  Z  value  for  ou^t  purposes 

-  Maximum  grid  dimension 

-  Minifflum  grid  dimension 

-  Maximum  R  value 

-  Amplilude  threshold  when  spectral  term  is  added 

-  Minimum  number  of  spectral  terms  in  a  cakualtkm 

-  Maximum  number  of  spectral  terms  in  the  calcualtkw 

-  Number  of  normal  tongiOklinal  output  points 

-  Number  of  detailed  longitudinal  output  poials 

-  Number  of  transverse  boundary  points  (induding  0) 

•  Total  number  of  output  points 

-  Maximum  transverse  index  for  output 

-  Maximum  array  index 

-  Maximum  transverse  army  index 
-IMAX-1 

-IMAX-2 

-  Starting  transverse  array  index 

-  Maximum  transverse  army  index 

-  Flag  for  restart  after  a  oontroUed  exit 

-  Counter  for  on-axis  output 

-  Counter  for  off-axis  0U90t 

-  Local  counter  for  fixing  output  points 

-  The  current  analysis  case 

-  Step  number 
-Pi 

-  2.0  •  Pi 

•  Initial  On-axis  pressure  (kPa) 

-  Angular  fundamental  frequency 

-  Nonlinearlity  parameter  (fi/A) 

-  1  +  BoverA/2 

•  Number  of  uniformely  spaced  radial  poiMs 

-  Radial  step  size  in  meters 


DOUBLE  PRECISION  K1 ,  K2,  K3,  KN1 ,  KN2,  KN3,  XN4 
DOUBLE  PRECISION  ZMIN,  ZMAX,  Threshold,  Pi,  TsraPi 
DOUBLE  PRECISION  PO,  Onega,  BoverA,  BoverZA,  Rnax,  Rgnirt,  Rgnax 
DOUBLE  PRECISION  OldOeltaR,  NewlOeltaR 

INTEGER  CASE,  MAXM,  MINM,  RESTART,  lOUTMAX,  NPOINTS,  NDPOINTS,  Nunber  BP 
INTEGER  TPOINTS,  COUNTER,  IMAX,  I,  NSTEP,  K,  L,  M,  IHAXMIN 
INTEGER  IHAXMAX,  IMAX4,  IMAXl,  Nuaber  RP 
LOGICAL  NEUffEST 

This  is  the  first  set  of  COMMON  blocks 


COMMON  /ALL/  NEWTEST 

COMKM  /CONSTANTS/  A,  0,  Freq,  CO,  TwoPi,  Pi 
COWON  /KVALUE/  K1 ,  K2,  K3,  KN1 ,  KN2,  KN3,  KN4 

This  common  is  used  for  indicces 


COMMON  /INDEXES/  IMAX,  M,  IMAXl 
Ihb  common  is  used  by  NBCT  and  the  main  program 


COWiON  /NEXTC/  OeltaR,  OeltaZ,  Rayleigh,  Ld,  AlphaD,  (Max 

At  tins  point,  take  the  maximum  and  minimum  Z  values  and  the  specified 
detailed  output  points  and  determine  the  points  to  ouqtut  data. 
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The  logic  is  simple.  Fm  calculate  the  location  of  eack  Bonnal  output 
poiaL  Then  make  a  loop  and  use  tke  loop  to  populate  the  ZOUT  anay.  As 
the  anay  is  populated,  check  to  see  if  a  detailed  output  point  has  been 
reached.  If  so,  add  it  to  the  amy 

Before  reacfing  the  input  data  write  one  as  the  restart  flag  so  that  the 
restart  option  can  be  eaercised.  if  desired. 

0PEN(S,F1LE  -  ’RESTART*) 

REAOfS,*)  Restart 
CL0SE(5) 

Set  the  matimuffl  array  index 
HARRAYINOEX  >  4000 

If  restarting  execution,  reset  the  current  parameter  vahies  to  the  ones  pment 
when  the  ran  was  tenninated. 


IF  (RESTART. EQ.1)  THEN 


299 

298 


If  restarting,  read  the  relevant  paramters. 
0PEN(5,FILE»’TEMP1 .DAT') 


REA0(5,*) 
REA0(S,*) 
REA0(5, 
READ(5, 
READ(5, 
REA0(5. 
READ(S, 
REA0(5, 
REA0(5, 
REA0(5, 
REA0(5, 
READ(S, 
REA0(S, 
REA0(S, 
READ(S, 
REA0(S, 
REA0(5, 
REA0<5, 
READ(5, 
REA0(5, 
READ(5, 
REA0(5, 
REA0(5, 
READ(5, 
READ(5, 
REA0(5, 
READ(5, 
REA0(5, 
REA0(5, 
REA0(5, 
REA0(5, 
READ(S, 
READ(5, 
READ(5, 
REA0(5, 
READfS, 
REA0(5, 
REA0(5, 
DO  299 


Case 

Threshold 
loutHax 
Rflsax 
Rgaln 
Raax 
AlphaO 
Alpha 
AiphaO 
PO 
RoO 
BoverA 
Bover2A 
0 
A 

DeltaR 
DeltaZ 
Zain 
Zaax 
Z 
K 
L 

NSTEP 
IMAXHAX 
IHAXMIN 
IMAX 
IIMX4 
K1 
K2 
K3 
KN1 
KN2 
KN3 
KN4 

Npolnts 
NDpolnts 
Nuaber  BP 
Tpoints 

9  299  Z  -  1.  NDpolnts 
REA0(5,*)  hi(I) 

CONTINUE 

00  298  1  >  0,  Nuaber  BP  -  1 
REA0(S,*)  8oundaryP(0,I),BoundaryP(1 ,1) 


CONTINUE 
READ(5,*) 
REA0(5,*) 


MlnH 

M 

MaxM 


READfS,*) 
REA0(5,*)  CO 
REA0(5,*)  Freq 


C4 


I* 

•tgh 
ReA0(5,*)  TwoPi 

Now,  read  (he  iniliaUy  caicualted  ouipuc  array. 


REAO(5.*)  Oae 
REAO(5.*)  Ray 
READ(5,»)  Ld 
READ(5  *)  PI 


DO  606  I  -  1 ,  Tpolnts 
«EAD(5,»)  Za/T(I) 

606  CONTINUE 
CL0SE(5) 

*  Next,  read  the  amyi  naed  in  the  cakulaiioa  G  and  H. 

CALL  REAOTEMP(G,  H,  MAXH,  IMAXMAX) 

■  Now,  open  the  old  output  filet  from  the  lermiaated  run.  set  the  poiate* 

*  to  the  end  of  the  file,  and  tetm  to  the  program. 

334  0PEN(11,FILE>'AHP2’  STATUS>’0L0') 

READO1,*,EN0>335)  OUMIY 

00  TO  334 

335  0PEW<12,PILE«’PHASEZ',STATUS*’0L0*) 
REA0f12.*,EN0«336)  OUIHY 

00  TO  335 

336  0PEM(13,F1LE=‘A1IPZR’  , STATUS* ’OLD* ) 

REA0(13  *,EN0*337)  OUNVy 

GO  TO  336 

337  0PEN(14,FILE*’PHASEZR'  .STATUS* ’OLD’  ) 
READ(14,*.EN0=338)  OUMIY 

GO  TO  337 

338  BACKSPACE(II) 

BACKSPACE(12) 

BACKSPACE!  13) 

BACKSPACE04i 
OUMIY  «  OUMIY  •  1  .0 

* 

*  Jump  to  program  exeoution. 

GOTO  333 
ENOIF 

• 

*  This  is  not  a  restart,  therefore  open  the  menu  program.  But  first  set 

*  restart  to  1  for  future  reference. 


OPEN  (5, FILE  »  ’RESTART’) 
Restart  ■  1 
WRITE(5,397)  Restart 
397  F0RIIAT(f4) 

CL0SE(5) 

*  Execute  the  mean  program. 


CALL  SYSTEU(’llenu.exe’) 


This  is  the  primary  code  input  sectioa.  For  the  integrated  Pascal  A  FORTRAN 
code,  the  file  INPUTDAT  is  the  transfer  point  between  the  two  codes. 


0PEN(5,FILE  »  ’INPUT. OAT’) 


REAOfS, 

READfS, 

REAOfS, 

REAOfS, 

READfS, 

REAOfS, 

REAOfS, 

REAOfS, 

READfS, 

REAOfS, 

READfS, 


Case 
Threshold 
loutMax 
Rgaax 
AlphaO 
PO 
RoO 

BoverA 

0 

A 

OeltaR 
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399 


REA0(5,«)  Ziin 

- 

Npoints 
NOpoints 

Do  399’ I  ■  1,  NOpoints 
REA0(5.*)  h)(I) 

CONTINUE 
NEAD(5,*)  MinM 
REAOis,*)  MaxM 
R£AD(5,»)  CO 
READ(5,*)  Freq 
REAIJ(5,»)  l^niber_BP 
Do  297  I  V  d  Numter_BP  -  1 
REAtHV)  BoundsiyP(04).Bound»fyP(l.I) 
297  CONTINUE 


REAOfS,*)  2m 
REA0(5,»)  Np 
REAOls,* 


CIOM  Che  file  and  calculate  some  of  the  paraaeten 


Ct.OSE(S) 

Set  the  boundary  regioa  Rgmin  to  10%  greater  than  the  physical  edge 
of  the  traiuducer. 


Rgaln  •  1 .100 
Calculate  ^2A  usiitg  B/A 


Bover2A  ^  1 .ooo  *  (BoverA/2.000) 
Set  the  values  of  Pi  and  2  Pi  for  the  program 


Pi  «  2.000  •  OASIN(I.OOO) 

TeroPi  *  2.000  ♦  Pi 

Set  Omega  and  the  attenuatioa  terms 


Oaega  «  TwoPi  •  Freq 

Alpha  •  AlphaO  ♦  (Fr««i»»2) 

AlphaO  •  Alpha  *  0 

Set  the  total  number  of  points 


TPOINTS  =  NPOINTS  ♦  NOPOINTS 
Now  open  the  output  files 


0PEN(11,Flle  «  ’AMPZ’l 
0PEN(12,Flle  »  ’PHASEZ’) 

0PEN(13,FILE  >  ’AMPZR’) 

0PEN(14,F1LE  «  *PHAS^') 

*  To  identify  the  output  later  include  the  input  data  a  file  header. 

00  237  I  «  11,14 

CALL  HEA0ER(1, Case, 0, A, PO,DeltaR,Zain,Zaax,Ioutliax, Npoints, NOpoints, Nuaber  BP.Tpolnts, 
ZU, MinM, UaxM, Freq ,BoundaryP) 

237  CONTINUE 

*  Set  Z  to  ^IN.  COUNTER  to  I,  and  populate  the  output  point  array. 

* 

DUHUY  >  (ZMAX  •  ZMIN) /NPOINTS 
Z  >  ZMIN 
COUNTER  -  1 
00  10  I  «  1,  TPOINTS 
Z  >  Z  ♦  OUMUY 

IF^((Z^GE.ZUi(^NTER))  .AND.  (COUNTER. LE. NOPOINTS))  THEN 

ZOUT(I)  >  ZU(COUNTER) 

COUNTER  COUNTER  ♦  1 
ELSE 


C-6 


20UT(I)  >  Z 
ENOIF 

10  CONTINUE 

Now  initialize  a  number  of  program  parameters. 

Fnt.  calcualte  the  K  values  used  in  the  <£flereacing  scheme. 


K1  -  -1.000/12.000 
K2  >  16.000/12.000 
K3  -  -30.000/12.000 
KN1  >  1.000/90.000 

KN2  >  -1.000  *  ((1.000/12.000)  ♦  (2.000/30.000)) 
KN3  «  (16.000/12.000)  ♦  (15.000/90.000) 

KN4  >  -1.000  *  ((30.600/12.000)  ♦  (2.060/9.000)) 

Then  the  rest 


lUMOlAX  -  XNT((Rgiax/OeltaR)  ♦  1.00-5) 

This  keep*  the  data  entiy  enon  from  btoving  iqt  the  code. 

IP  (IIMXUAX.GT.MARRAYINOEX)  THEN 
IHMCIIAX  >  MMWAYINOEX 
ENOIF 

Continue  srith  the  calculalkM. 


IMAXMIN 

>  INT((Rgain/OeltaR)  * 

1 .Od-5) 

IMAX 

=  IMAXMIN 

IMAX4 

-  IMAX  -  4 

K 

3  1 

L 

3  1 

NSTEP 

3  0 

Rayleigh 

3  Onega  *  A  *  A/ (2. 000 

•  CO) 

For  later  use  save  the  initial  vakie  of  OehaR 


OldOeltaR  >  OeltaR 

OeltaR  >  OeltaR  *  D5QRT(Rayleigh/D) 

Raax  -  OSQRT(Rayleio1)/0) 

Use  0.4  in  the  DeltaZ  oonvenion  process  to  make  the  calculatkxB  more  suble. 

DeltaZ  «  0.4d0  *  (OeltaR  **  2) 

Z  «  -DeltaZ 

Calculate  the  initial  thaoontinuity  distance. 

Ld  ■  RoO  •  (C0*»3)/(Bover2A  •  Onega  •  PO) 

H  =  Minn 

Now  run  the  initializntioo  programs.  Remember  the  initializadon  bat  to 
be  at  the  matimnm  grid  points.  Fast  sei  the  transverse  boundary  oondiiioa. 

CALL  SET_BOUNIlARY(A,  Number_BP.  Radial_P.  BoundatyP.  NewDellaR,  Number_RP,  OldDeluR,  IMAXMAX) 


CALL  IN1TIAL(G,  H,  IMAXMAX,  MAXM,  OeltaR) 

*  This  is  the  rentry  point  for  a  restart. 

333  CONTINUE 

IMAXI  *  INAX  -  1 

*  To  allow  program  restart,  set  the  ooatrol  break  option.  Once  Ctrl  C  is 
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prcMed,  the  prognm  will  jump  to  the  desipiated  point,  once  detected. 


NEWTEST  «  .FALSE. 
OPTION  BREAK(NEWTEST) 

Set  the  m^or  prognm  loop 


100  NSTEP  >  NSTEP  *  1 
2  >  Z  4-  OeltaZ 
IF  (Z.6T.20UT(K))  THEN 
2  -  2DUT(K) 

ENOIF 

CALL  NEXT(G,H) 


Ihu  simple  four  Gac  section  of  code  floats  the  traasvene  boundary. 
Chlculate  the  ampftude  four  cells  from  the  rigid  boundary.  If  it  has 
emeeded  the  threshold,  eipand  if  necessary.  This  worlct  because  the  lowest 
index  value  oorreapomls  to  the  lowest  frequency. 


IF  (IMAX.LT.IMAXMAX)  THEN 

• ' )  >^ ( 1 1 ) )  .GT .Threshold )  THEN 

IF  (IHAX.GE.IMAXIIAX)  IMAX  »  IMAXliAX 
IIMXI  >  IMAX  •  1 
IMAX4  >  IMAX  -  4 
ENOIF 
ENOIF 

Now  see  whether  the  on-axis  output  is  ttecesssary  by  checking  to  see  if  the 
value  corresponds  to  an  element  of  zout(k). 


IF  (Z.EQ.ZOUT(K))  THEN 


At  each  output  point,  write  the  on-axis  ampfitude  and  phase  spectra. 


500 


1009 


WRITEfll ,1009)  M 
WRITEM2,1009)  M  ,  2',  (PHASE(G<d,li) ;H( 
MUTEl*,^)  Nstep,  z,  k,  ■,  Ihax 
FORMATC  Nstep  «  ’,16,’  Z  »  ’,F10.6,’  k 
K  «  K  +  1 


’.13,’ 


O.M) 


ENOIF 

F0RMAT(1P,14,’  ’ ,500E16.4E4) 
IF(Z.EQ.ZU(L))  THEN 


’,13,’  max  »  ’,14) 


If  z  b  at  the  pont  z(1)  the  increment  the  index  in  the  zu(l)  array. 


IF(L.LT.NOPOINTS)  THEN 
L  »  L  ♦  1 
ELSE 

L  «  NOPOINTS 
ENOIF 


Thb  b  the  place  to  output  the  spectra  to  ampr  aud  phaser. 


2000 


00  2000  INEW  *  0,  lOUTMAX 
1IIIRITE(13,1009)  M,  2,  (OldOeltaR 
«IIRITE(14,1009)  M,  2,  (OldOeltaR 
CONTINUE 


INEW) ,  (AMPLTU(G(INEW,N)  ,H(1NEW,N) )  ,N^,H) 
INEW),  (PHASE(G(INEW,N),H(1NEW,N)),N-0,H) 


ENOIF 

IF  (M.LT.MAXM)  THEN 
IF^AMPI^TU(G(0,M),H(0,M)). 


GT. Threshold)  THEN 


ENOIF 

ENOIF 

IF(K.EQ.(TP0INTS+1))  GO  TO  115 


*  If  has  been  pressed  the  go  to  90S. 

9 

IF(NEWTEST)  GOTO  905 
GOTO  100 
115  CONTINUE 


If  Soidied.  tbeo  react  resurt  to  0. 


OI»OI(5.FILE*’RESTART‘ ) 
RESTART  >  0 
«(RITE(5,«)  RESTART 
CL06E(5) 

Now  jump  to  906. 


GOTO  906 
905  CONTINUE 

*  If  reataitiog,  then  write  the  relevant  paranten. 


199 


296 


h(5/I 

WRITEfS, 

WRITEfS, 

WlTEfS. 

NRlTEfS, 

WUTE(5, 

MRITEfS, 

MUTE(5, 

MRITEfS, 

MRITEfS, 

MRITEfS,* 
MRITEfS,*)  RoO 
MRITEfS,*)  BoverA 
MRITEfS,*)  Bover2A 
MRITEfS,*)  0 
MRITEfS,*)  A 
MRITEfS,*)  OeltnR 
MRITEfS,*)  OeltaZ 
MRITEfS,*)  Znin 
MRITEfS,*)  Znax 


Cnsn 

Threshold 

loutHax 

Roanx 

Rgaln 

rSw 

AlphnO 
Alpha 
AlphaO 
PO 


MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
WIITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*)  K3 
MRITEfS,*)  KN1 
MRITEfS,*)  KN2 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 
MRITEfS,*) 


Z 
K 
U 

MSTEP 

INAXMAX 

IHAXMIN 

INAX 

IMAX4 

K1 

K2 


KN3 

KN4 

NMlntS 
NOpolnts 
Ntwber  BP 


- 

MRITEfS,*)  Tpoints 
00  199  I  ■  1,  NOpolnts 
MRITEfS,*)  ZUfI) 

CONTINUE 

00  296  I  >  0,  NtMiber  BP  •  1 
MRITEfS,*)  BoundaryPfO,!),’ 
CONTINUE 
MRITEfS,*)  MlnM 
MRITEfS,*)  H 
MRITEfS,*)  MaxH 
MRITEfS,*)  CO 
MRITEfS,*)  Freq 
MRITEfS,*)  Oeeaa 
MRITEfS,*)  Rayleigh 
MRITEfS,*)  Ld 


MRITEfS,*) 

MRITEfS,*) 


Pi 

TwoPl 


’ ,BoundaryPf1 ,1) 
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Now  write  the  ouqait  array  that  was  iniiialh  calculated. 


00  406  I  ■  t ,  Tpoints 
lll(RITE(5,»)  ZOUT(I) 

406  CONTINUE 
CLOSE (5) 

*  Next  write  the  arrayi  used  for  G  and  H  in  the  calculatioa 
CALL  INRITETENP(6,  H,  liAXM,  IMAXUAX) 

906  CONTINUE 

miRITEf*,*)  ’PROGRAM  EXECUTION  COtTLETEO* 

*  Now  dose  the  output  Ska 

CL08E(11) 

CL0SE(12) 

aOSEfIS) 

CL0SEO4) 

STOP 

END 
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SUBROUTINE  HEADER ( N ,  Case , 0 ,  A , PO , Del t aR , ZMln , Zma» ,  ToutHax , Npoints , NDpolnts , Nuaber  BP , Tpoints , 
ZU ,  MtnH ,  tIaxM ,  Freq ,  BoundaryP) 

*  This  subroutine  writes  a  header  on  each  output  file 

* 

DOUBLE  PRECISION  Zaln,  Zaax,  DeltaR.  0,  A,  Freq,  ZU(50).  PO.  BoundaryP(0:1 ,0:100) 

INTEGER  CASE,  I,  Npoints,  NOpoints,  Tpoints,  liinll,  Maxit,  N,  loutUax,  Nuaber  BP 
LOGICAL  NEWTEST 
COMMON  /ALL/  NEWTEST 

• 

*  This  is  the  ooUectioo  of  write  statemeats 


1002 


1003 

1000 


WRITEfN,*)  CASE 
WRITE(N,1000)  0 
NRITE(N,1000)  A 
WRITE(N,1000)  PO 
WRITE(N,1000)  OeltaR 
WRITE(N,1000)  Freq 
\NRITE(N,1000)  Zain 
vmiTE(N,1000)  Zaax 
miTEfN,*)  loutMax 


WRITE(N,» 

WRITE(N,* 

WRITE(N,» 

NRITE(N, 


Npoints 
NOpoints 
Tpoints 
MlnM 
miTElN,*)  MaxM 
DO  1002  1  =  1, NOpoints 
WniTEfN.IOOO)  ^(I) 

CONTINUE 

DO  1003  I  >  0,  Nuaber  BP  •  1 
WRITF(N,»)  BoundaryP(0,I),* 
CONTINUE 

F0RMAT(1P,E15.6c3) 


RETURN 

END 


' ,BoundaryP(1 ,1) 
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SUBROUTINE  READTEUP(G,H.MAXM.IMAX) 

DOUBLE  PRECISION  G(0:4060.0:30) ,  H(0:4000,0:30) 
INTEGER  MAXV.IMAX 

This  subroutine  reads  fUes  TempZdai  through  Temp54aL 


608 

607 


308 

307 


6t8 

617 


318 

317 


0PEN(5,FILE»’TEMP2.0AT' ) 

00  ear  i>o.maxh/2 
00  608  J  *  O.IUAX 
REA0(5,*)  G(J,I) 
CONTINUE 
CONTINUE 
CLOSE(S) 

OPEN  ( 5 ,  FILE=  ’  TEIIP3 .  OAT  ■ ) 
00  307  I-(IIAXM/2)4^1,IIAXH 
DO  308  J  *  O.IHAX 
READ(5,*)  6(J,I) 
CONTINUE 
CONTINUE 
CL0SE(5) 

0PEN(5,FILE=*TEMP4.DAT’) 
00  617  I«0,IIAXM/2 
00  618  J  >  O.IMAX 
READ(5,*) 

CONTINUE 

CONTINUE 

CLOSE(S) 

OPEN ( 5 , F I LE= ■ TEMPS . DAT ■ ) 
00  317  Is(MAXM/2)-»1,UAXM 
00  318  J  >  O.IMAX 
READ(5,»)  H(J,I) 
CONTINUE 
CONTINUE 


CL0SE{5) 

RETURN 

END 
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SUBROUTINE  «IRITETElP(G,H,IMXH,IltAX) 

DOUBLE  PRECISION  G(0:4006,0:30) ,  H(0:4000 
INTEGER  UAXM.IMAX 

This  subroutine  writes  Teinp2xiat  through  TempSiial 


SO 


90 

100 


190 

200 


290 

300 


390 

400 


F0RIMT(1P,4030.15) 

OPEN ( 5, F I LE>’ TE1IP2.DAT' ) 
00  100  I»0,IMXM/2 
00  90  J  *  O.IMAX 
1WITE(5,‘)  0(J,I) 
CONTINUE 
CONTINUE 
aOSE(5) 

0PEN(S,FILE-‘TEMP3.0AT' ) 
00  200  I«(IIAXM/2)*1,IMXM 
DO  190  J  «  O.IMAX 
INRITEtS,*)  G(J,I) 
CONTINUE 
CONTINUE 
CLOSEtS) 

0PEN(S,FILE>>TEMP4.0AT‘) 
00  300  lsO,ll^Xll/^ 

DO  290  J  *  O.IHAX 
WRITE(5,»)  H(J,1) 
CONTINUE 
CONTINUE 
CL0SE(5) 

0PEN(S,FILE»’TEMP5.0AT' ) 
00  400  I>:(IU0(M/2)41,IIAXH 
DO  390  J  >  O.IMAX 
»KRITE(5,»)  H(J,1) 
CONTINUE 
CONTINUE 
CL0SE(5) 

RETURN 

ENO 


.0:30) 
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100 

110 


SUBROUTINE  SET  BOUNOARY(A.  Nuaber  Bf>,  Radial  P,  BoundaryP,  NewDeltaR,  Nuaber  RP, 

oidOaltaR,  laax)  ~ 

This  subroutine  cilcuiates  the  initial  nonnalized  picssuic  distrbutioa 
Declared  variables 


A 

Nuniber_BP 

Number_RP 

BoundaryP 

OldDeltaR 

NewDdtaR 

Radius 

Radial_P 

Localindei 

I 

Imaz 

m 

b 


s  Transducer  radhis 
»  Number  of  eipcrimental  values 

>  Number  of  unifarm^  spaced  radial  points 
m  Eiperimental  point  2-D  amy 

a  Normalized  step  size 
K  Radial  step  size  (meters) 

>  Radial  location  (meters) 

•i  Uniform^  spaced  normalized  pressure  values 

-  Array  in^i  for  incrementing  BoundaryP 
m  Loop  indez 

»  Maximum  number  of  transverse  terms 

-  Sfope 

«  Y  intercept 


INTEGER  Nuaber  BP,  Localindex,  I,  Nuaber  RP,  laax 

DOUBLE  PRECISION  Radial  P(0:4000),  BoundaryP(0:1 ,0: 100) ,  NewDeltaR 

DOUBLE  PRECISION  OldDelTaR,  Radius,  a,  b,  a 

Set  number  of  radial  points.  This  value  is  passed  to  main  program. 


Nuaber  RP  -  IDINT(1 .OdO/OldOeltaR) 

IF  (DMDD(1  .OdO, OldDeltaR). NE.O.OdO)  Nuaber_RP  =  Nuaber_RP  *  1 

Given  the  number  of  radial  points,  set  the  pressure  values.  For  a 
uniform  distribution  set  the  pressure  to  1  and  go  to  the  end. 


IF  (Nuaber  BP.EQ.1)  THEN 
DO  100  I*a  0,  Nuaber  RP 
Radial  P(l)  >  I.OdO 
CONTINUE" 

DO  110  I  *  Nuaber  RP  'fl ,  laax 
Radial  P(l)  >  OTOdO 
CONTINUE" 

GOTO  500 
ENDIF 


Now  if  we  don't  have  a  uniform  distrlwtioo,  construct  a  profile 
based  on  the  eiperimental  points.  Set  the  on  axis  point  (0), 
the  local  index,  and  tadial  step  size  (in  meters). 


Radial  P(0)  >  BoundaryP(1 ,0) 

Localindex  >  1 

NewOeltaR  »  OldDeltaR  *  A 


Now  build  the  amy 


DO  200  I  >  1,  Nuaber  RP 
Radius  «  1  *  NewOeltaR 


Check  whether  the  radius  is  <  Boundar)'P(0JocaEndex).  If  so, 
interpolate  to  rletcrmine  Radial_P(I). 


210  CONTINUE 

IF  (Radius. LE.BoundaryPfO, Localindex))  THEN 
a  *  (BoundaryP(1 (Localindex)  •  BoundaryP(1 .Localindex  -  1)) 
a  =  a/ (BoundaryP(0, Localindex)  -  BoundaryP (0, Local Index  •  1)) 
b  *  BoundaryP(1 (Localindex)  -  (a  *  BoundaryPiO, Localindex)) 
Radial_P(I)  ^  (a  *  Radius)  *  b 


If  not,  increment  the  index  and  loop  back. 
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IF  (Locallndex  .LT.  (Njaber  BP  •  1))  THEN 
Locallndex  -  Localindex  ■»~i 
GOTO  210 
ELSE 

Radlal_P(I)  ■  BoundaryP(1 (Locallndex) 
ENOIF 
ENOIF 

200  CONTINUE 

DO  220  I  *  Nuaber  RP  >1 ,  laax 
Radial  P(I)  >  OTOdO 
220  C0NTINUE~ 

500  CONTINUE 

• 

*  Now  write  the  results  to  check  the  prooes. 

• 

0PEN(5.File  •  ’TeapFile.dat’ ) 

Wrlta(S,*)  Nuaber  RP 
Writers,*)  NewOelfaR 
Nrlta(5,*)  A 
Writers,*)  Nuaber  BP 
Writers,*)  OldOelTaR 
Do  300  1  ■  0,  laax 

RADIUS  >  I  *  NewOeltaR 
writers,*)  Radius,'  Radial  P(l) 

300  continue  ■ 

CLOSE(S) 

*  Slop  and  return. 


RETURN 

END 
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SUBROUTINE  IN1TIAL(Q,  H,  IMAX.  MAXM,  OeltaR) 

This  subroutine  Initializes  the  coeffickats  G  and  H  at  the  transvcise  bouodaty. 


DOUBLE  PRECISION  0(0:4000,0:30),  H(0:4000,0:30) ,  U,  OeltaR 
DOUBLE  PRECISION  INIAMP,  INIPHS,  A,  0,  Fre<|,  CO,  TwoPl,  PI 
INTEGER  UAXH.  IMAX,  I,  N 
LOGICAL  NEWTEST 
COIHON  /ALL/  NEWTEST 

COMKM  /CONSTANTS/  A,  0,  Freq,  CO,  TwoPl,  PI 
Since  only  a  single  hamonic  eiists  at  the  start,  initaEze  C  and  H  of  (Ll). 


10 


00  10  I>0,IMAX 
U  -  1  *  OeltaR 

IF  (Radial  P(I).NE. 0.000)  THEN 
0(1,1)  «~Radial  P(1)*OCOS(INIPHS 
h(1,1)  -  Radlal~P(I)*OSlN(lNIPHS 
ELSE 

G(I,1)  >  0.000 
h(I,1)  >  0.000 
BIOIF 
CONTINUE 


For  safety’s  sake  set  all  other  teims  to  0 


00  20  1  >  0,  IMAX 
00  30  N  =  2,  IMXM 
G(I,N)  «  6.000 
H(I,N)  -  0.000 
30  CONTINUE 
20  CONTINUE 
RETURN 
ENO 
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D0U8I.E  PRECISION  FUNCTION  INIPHS(U) 

Hiit  procedure  ad^ts  tbe  phac  of  the  input  pressure  (root  so  that  it 
appears  to  be  oomming  from  a  curved  surfKe. 


DOUBLE  PRECISION  A,  D,  Rvalue,  U,  Freq,  CO,  TwoPi,  Pi,  K 
DOUBLE  PRECISION  OeltaR,  OeltaZ,  Rayleigh,  Ld.  AL^aO  Raax 
LOGICAL  NEWTEST 
COMMON  /ALL/  NEWTEST 


COMMON  /CONSTANTS/  A,  0.  Free,  CO,  T««Pi,  PI 

COMMON  /NEXTC/  OeltaR,  Deltai,  Rayleigh,  Ld,  AlphaO,  Raax 


Convert  from  normalized  radiM  ooordinaies  to  real  coordinates. 


Rvalue  >  A  *  U  *  DSQRT(0/(4.000  *  Rayleigh)) 
Calciilate  the  wave  number  baaed  oe  the  input  frerpieacy 

K  ■  TwoPl  •  Freq/CO 
Phase  Factor 

INIPHS  >  K  *  (0SQRT(0*0  *  Rvalue* Rvalue)  -  0) 

return 

END 
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SUBROUTINE  NEXT(G,H) 

This  subroutine  calculates  G  and  H  for  the  next  longiiudiaai  distance  step. 


DOUBLE  PRECISION  0(0:4000,0:30),  H(0:4000,0:30) 

DOUBLE  PRECISION  GNEW(0:4000) ,  HNEW(0:4000),  NORM 

DOUBLE  PRECISION  GFN(0:4000),  HFN(0:4000).  00LO(O:4O0O) ,  Raax 

DOUBLE  PRECISION  SUMO,  SUMH,  Eps,  Rayleigh,  A,  0,  Freq,  CO,  03 

DOUBLE  PRECISION  01,  02,  OeitaR,  DeltaZ,  Ld,  AlphaD,  TwoPl,  Pi 

INTEGER  M,  IMAX,  LMAX,  IMAX1 

LOGIOAL  N^EST 

OOlllON  /ALL/  NEWTEST 

COMMON  /CONSTANTS/  A,  0,  Freq,  CO,  TWoPi,  Pi 

COMMON  /NEXTC/  DeltaR,  Deltai,  Rayleigh,  Ld,  AlphaO,  Raax 

COIMON  /INDEXES/  IIMX,  M,  1MAX1 

*  Set  the  local  conveigence  oonslants 

Eps  >  1.0d-5 
LMAX  -  20 

*  These  constants  are  calculated  outside  the  loop  to  speed  np  the  prograna. 

* 

DO  60  N=1,M 

01  >  DeltaZ  *  N  *  0/(2.0D0  *  Ld) 

02  «  1.000  /(1 .000  *  (DeltaZ  *  AlphaO  *  N  •  N)) 

03  3  DeltaZ/ (4.000  *  N  *  OeitaR  *  OeitaR) 

*  This  loop  calculates  the  non-Enear  term  oontrSMition. 

• 

DO  10  I=0,IMAX1 
flOLO(I)  =  0(1, N) 

GFN(I)  -  G0L0(I)«C1‘SUMG(G,H,N,M,I) 

HNEW(i)  >  H(I,N) 

HFN(I)  -  HNEW(I)-fC1*SUMH(0,H,N,H,I) 

10  CONTINUE 

• 

*  L  is  the  convergence  counter.  If  L  exceeds  20.  the  ptognm  aborts. 

L  »  1 

*  Apply  the  seven  point  Lapacian 

20  CALL  C0NV0L(GNEW,GFN,C2,C3,HNEW) 

CALL  C0NV0L(HNEW,HFN,C2,-C3,QNEW) 

*  If  the  step  to  step  error  exceeds  the  oonvergenoe  bctor.  continue  to  iterate. 

* 

DO  50  I-0,IMAX1 

NORM  -  6aBS(GNEW(I)  •  OOLD(I)) 

IVAL  *  I 

IF  (NORM.GT.EPS)  GOTO  77 
50  CONTINUE 

• 

*  If  the  setp  to  step  error  is  allowable,  jump  to  the  next  step. 

• 

GO  TO  54 
77  CONTINUE 

DO  59  I>0,IMAX1 
GOLO(I)  =  GNEW(I) 

59  CONTINUE 

• 

*  If  L  >  20,  then  terminate  the  program  and  write  the  diagnostics. 


IF  (L.GT.LMAX)  THEN 

NRITE(»,»)  'ITERATION  FOR  N  =  ’,N,'  FAILED.  PROGRAM  ABORT’ 
WRITE(*,885)  NORM,  Eps,  L,  IVAL 
885  F0niAT(1p,2d15.6,h5) 

STOP 

ENDIF 
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Increment  L  and  continue  to  iterate 


L  =  L+1 
GO  TO  20 

Set  G  and  H  to  the  new  values 


54 


00  53  1  «  O.IUAXI 

Q(I,N)  >  GNEW(I) 
»  HNEW(I) 

CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 
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SUBROUTINE  C0NV0L(Z,X,C2,C3,Y) 

This  subfouunc  perfonns  the  seven  point  lapUcian  operator  at  each  point  in  the  Anay  Y. 


DOUBLE  PRECISION  Z(0:4000),  Y(0:4000),  X(0:4000),  C2,  C3,  K1,  K2,  K3 

DOUBLE  PRECISION  KN1 ,  KN2,  KN3,  KN4 

INTEGER  IMAX,  H,  1MAX1 

LOGICAL  NEWTEST 

COtMON  /ALL/  NEINTEST 

COMMON  /INDEXES/  IMAX.  M,  IUAX1 

COMMON  /KVALUE/  K1 ,  K2,  K3,  KN1  ,  KN2,  KN3,  KN4 

Z  «  The  output  value  which  ia  G  or  H 

X  s  The  previous  giidpoint  value  which  it  G  or  H 

Y  -  The  current  grid  value  which  is  H  or  G 

C2  >  IjODO /(IjODO  +  (DeltaZ  *  Ai^haO  •  N  *  N)) 
a  -  D  *  OeltaZ/(4.0D0  •  ZD  «  N  *  DeltaR  *  DeltaR) 


Calcualte  the  oa-tak  term,  and  the  next  two  tenns  moving  away  frcm  the  aiia 


Z(0)  »  C2  *(X(0)+  C3  •  (KN1  •  (Y(3)  ♦  Y(3)) 

♦  KN4  •  Y(0))) 

Z(1)  »  C2  •(X(1)+  63  *  (KN1  •  (Y{4)  ♦  Y(2)) 

♦  «M  *  Y(1))) 

Z(2)  =  C2  •(X(2)+  C3  •  (KN1  •  (Y(5)  ♦  Y(1)) 

♦  KN4  •  Y(2))) 

Next  compute  the  terms  up  to  the  IMAX  -  4 


KN2  •  (Y(2)  ♦  Y(2))  ♦  KN3  •  (Y(1)  ♦  Y(1)) 
KM2  •  (Y(3)  ♦  Y(1))  ♦  KN3  •  (Y(2)  +  Y(0)) 
KN2  •  (Y(4)  ♦  Y(0))  ♦  KN3  *  (Y(3)  ♦  Y(l)) 


DO  10  I  =  3,  IMAX  -  4 

2(1)  »  C2  C3  •  (KN1  •  (Y(I+3)  ♦  Y(I-3))  ♦  KN2  •  (Y(I+2)  *■  Y(I-2)) 

+  KN3  •  (Y(l+1)  ♦  Y(I-1))  +  KN4  •  Y(I))) 

CONTINUE 


End  by  calcuailing  the  terms  at  the  upper  boundary 


Z(IMAX-3)  =  C2  «  (X(IMAX-3)* 
+  Y(IlM0(-4))  ♦  K3  • 
Z(IUAX-2)  *  C2  •  (X(IMAX-2) 
Z(IMAX-I)  >  C2  •  (X(IMAX-I) 


C3  •  (K1  •  (Y(IMAX-I)  ♦  Y(IMAX-5))  ♦  K2  •  (Y(IMAX-2) 
Y(IMAX-3))) 

♦  C3  •  (Y(IMAX-I)  •  2.000  •  Y(lMAX-2)  ♦  Y(lMAX-3))) 

♦  C3  •  (Y(lMAX-2)  -  2.000  •  Y|1MAX-1))) 


Return  to  the  mam  program 


RETURN 

ENO 
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Hiis  pair  of  fiiiictioas  is  used  to  calculate  linearized  sunt  teno  at 
each  longitudiiial  step 


10 


20 


DOUBLE  PRECISION  FUNCTION  SUIIG(G,H,N,M,I ) 
This  functioo  calcualtes  the  lineaiized  contriubtioa  for  G 


DOUBLE  PRECISION  6(0:4000,0:30),  H(0:4000,0:30) ,  SUM 

INTEGER  K.  N,  M,  I 

LOGICAL  NEWTEST 

COMMON  /ALL/  NEWTEST 

SUM  «  0.000 

00  10  K  *1,(N-1)/2 

SUM  »  SUM^(I,N-K)*Q(I,K)-H(I,N-K)*H(I,K) 

CONTINUE  /  t  .  » 

IF  (MOO(N,2).EQ.O)  THEN 
SUM  »  SUM  ♦  (0.5*(G(I,N/2)**2  -  H(I,N/2)**2)) 

ENOIF 

00  20  K=M,(N4-1),-1 

SUM  *  SUM-Q{I,K-M)*G(I,K)H(I,K-N)*H(I,K) 

CONTINUE 


SUMG  »  SUM 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  SUIW(6,M,N,M,I) 

*  This  function  calcualtes  the  linearized  contribution  for  H. 

• 

DOUBLE  PRECISION  G(0:4000,0:30) ,  H(0:4000,0:30) ,  SUM 

INTEGER  K,  N,  M,  I 

LOGICAL  NEVfTEST 

COMMON  /ALL/  NEWTEST 

SUM  »  0.000 

00  10  K>1,(N-1) 

SUM  *  SUM+G(I,N-K)*H(I,K) 

10  CONTINUE 

00  20  K=M,fN+1),-1 

SUM  *  SUM+G(I,K)*H(I,K-N)-H(I,K)*G{I,K*N) 

20  CONTINUE 

• 

SUWi  -  SUM 

RETURN 

END 
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The  following  funciioos  are  used  to  calculate  the  amlitude  and 
phaae  of  the  ipectfal  teims  givcii  the  local  vakies  of  G  and  11. 


DOUBLE  PRECISION  FUNCTION  AMPLTU(G.H) 

This  function  cakuahes  the  nonnalized  pnesure  amplitude  of  the  n*  banncnic. 

DOUBLE  PRECISION  G,  H 
LOGICAL  NEWTEST 
COMMON  /ALL/  NEWTEST 
AMPLTU  ■  0SQRT(6*G«H*H) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  PHASE(G,H) 


This  function  returns  the  nth  harmonic’s  phase  in  degrees. 


DOUBLE  PRECISION  G,  H,  MAG 
LOGICAL  NEWTEST 
COMMON  /ALL/  NEWTEST 

*  Calculate  the  magnitude  of  the  specual  lenn 

• 

MAG  »  OSQRT(G*G«R*H) 

IF  (MAG.EQ.O.)  THEN 
PHASE  »  0.000 
ELSE 

PHASE  =  0ATAN2(H,G)*90.000/0ASIN( 1.000) 

ENOIF 

« 

*  Make  sure  the  phase  is  a  number  between  0  and  3d0  degrees. 

• 

IF  (PHASE.GE.-ia0.00)  GO  TO  20 
10  PHASE  «  PHASE4-300.000 

IF  (PHASE. LE. -180.00)  GO  TO  10 
20  RETURN 

END 
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A  stroboscopic  Schlieren  image  of  ultrasonic  radiation  from  a  plate  immersed  in  water  is 
presented  and  is  compared  with  the  calculated  radiation  pattern  for  the  plate  in  which  a 
leaky  Lamb  mode  has  been  generated.  The  phase  relationship  of  the  radiation  lobes  created 
at  the  top,  bottom,  and  narrow  end  surfaces  of  the  plate  are  discussed. 

Keywords:  leaky  Lamb  modes;  Schlierert  imaging;  radiation  pattern  phase 


It  was  recently  shown'  ^  that  Lamb  waves  in  a  solid  plate 
immersed  in  water  radiate  some  of  their  energy  from  the 
end  of  the  plate.  This  radiation  from  the  end  of  the  plate 
consists  of  an  odd  number  of  lobes  when  the  plate  carries 
a  symmetric  Lamb  mode  and  an  even  number  of  lobes 
when  the  mode  is  antisymmetric.  The  previously  published 
Schlieren  images  produced  by  continuous  waves  and 
theoretical  calculations  of  the  relative  intensity  distribution 
between  the  lobes^,  as  well  as  hydrophone  measurements* 
do  not  indicate  the  respective  phase  of  the  radiation. 
Moreover,  the  use  of  continuous  waves  incident  from  the 
water  onto  the  plate  to  generate  a  given  leaky  Lamb 
mode  also  produces  specular  and  non-specular  reflections, 
first  observed  by  Schoch^,  and  later  explained  and 
described  by  Neubauer*,  Bertoni  and  Tamir’,  and 
others*'^.  When  these  reflections  of  the  incident  continuous 
wave  occur  near  the  end  of  the  plate,  they  may  interfere 
with  other  possible  reradiations  from  the  end  of  the  plate. 
The  Schlieren  images  published  by  Zhu  et  al.'  clearly 
show  reflections  from,  and  transmissions  through,  the 
plate  which  are  caused  by  such  interferences.  They  are 
not  specular  or  non-specular  reflections  since  the 
radiation  is  not  taking  place  at  the  expec'ed  Lamb 
angle. 

However,  the  numerical  approach  used  to  find  the 
intensity  distribution  of  the  radiation  at  the  end  of  the 
plate^  does  contain  information  from  which  relative 
phases  of  the  waves  in  the  lobes  can  be  determined. 
Furthermore,  the  numerical  method  can  be  used  to 
describe  radiation  from  the  top  and  the  bottom  cf  the 
plate,  near  its  end,  which  does  not  include  specular 
and  non-specular  reflections.  The  latter  are  simply  a 
consequence  of  a  particular  technique  of  exciting  Lamb 
modes  and,  thus,  may  be  considered  not  to  be  a  part  of 
the  radiation  from  a  Lamb  mode  per  se. 

This  paper  is  concerned  with  the  phase  and  the 
amplitude  distribution  of  the  radiation  occurring  at  the 


end  of  the  plate  as  well  as  from  the  top  and  the  bottom 
near  the  plate's  end.  The  interference  by  incident 
or  reflected  waves  has  been  eliminated  in  both  the 
mathematical  model  and  the  experimental  arrangement 
described. 


Calculational  approach 


The  basis  for  the  calculations  are  the  finite  different 
methods  discussed  by  Bond*  and  Marker’,  adopted  for 
the  present  case  of  energy  transfer  from  the  particle 
displacement  at  the  end  surfaces  of  an  aluminium  plate 
to  the  surrounding  water.  This  particle  displacement  can 
be  calculated  for  a  plate  vibrating  in  a  Lamb  mode  which 
is  excited  by  a  longitudinal  ultrasonic  wave  in  water, 
incident  at  the  appropriate  Lamb  angle. 

Assuming  that  the  long  dimension  of  the  plate  is  the 
x-direction,  the  thickness  of  the  plate  is  in  the  r-direction 
and  the  Lamb  wave  propagates  in  the  x-direction,  then 
the  wave  motion  in  the  homogeneous  solid  plate  can  be 
expressed  as 


{/.  +  p)W  u  -I-  = 


(1) 


Here  A  and  p  are  the  Lame  constants,  p  is  the  density  of 
the  plate  material,  and  u  is  the  displacement  vector  which 
has  an  x  and  y  component  and  is  time  dependent.  The 
>-direction  can  be  omitted  assuming  the  plate  to  be  wide 
enough  to  eliminate  any  y-dependence.  With  the  velocities 
ofthe  longitudinal  and  shear  waves  given,  respectively  by 


=  (/.  +  p),  p 


(2) 


=  p!p  (3) 

\ 

The  general  equations  for  the  particle  displacements  in 
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^  noo  c.' 


■1  n 


the  V  and  r  directions,  u  and  i,  rc>pccli\cl>.  are  then 

“ii  =  p'^xx  +  <r’  '  '•'•tj,.-  + 

f.l  =  P't:.-  +  (P*  -  <5) 

the  subscripts  indicate  derivatives. 

The  Lamb  wave  to  be  examined  is  produced  by 
allowing  an  ultrasonic  beam  of  frequency  /  (with  / 
between  2  and  10  MHz)  to  be  incident  at  the  appropriate 
angle  from  the  water  on  to  the  plate  so  that  the  desired 
Lamb  mode  is  set  up  in  the  plate.  Other  parameters 
needed  for  a  numerical  analysis  are  the  beam  profile, 
beam  width,  plate  thickness  <i,  and  distance  from  the 
point  of  beam  incidence  on  the  plate  to  the  end  of  the 
plate.  In  the  present  case  the  plate  was  0.9  mm  thick,  and 
the  distance  beam  to  the  end  of  plate  could  be  varied. 
The  beam  profile  was  assumed  to  be  Gaussian. 

Equations  (4)  and  (5)  tan  then  be  solved  by  the  FDM 
with  incremental  steps  of  .v  of  about  0.1/  and  the 
increments  of  t  are  selected  so  that 

At  ^  A.x(p^  4-  s^y  *  ^  (6) 

If  the  size  of  the  time  and  space  increments  are  selected 
so  that  Equation  ( 6)  is  satisfied,  the  accuracy  and  stability 
of  the  calculation  w  ill  be  assured 


Numerical  calculations 

The  above  method  was  used  to  calculate  the  radiation 
pattern  at  the  end  of  the  plate  vibrating  in  the  A, -Lamb 
mode.  The  frequency  of  the  continuous  wave  w  as  assumed 
to  be  3.29  MHz  and  the  plate  thickness  0.9  mm.  resulting 
in  an yif-value  of  2.96  MHz  mm.  The  A, -Lamb  mode  will 
be  excited  for  these  parameters  if  the  angle  of  incidence 
in  water  is  13.6  . 

This  set  of  parameters  is  used  as  the  input.  The  output 
is  an  array  of  calculated  sound-intensity  values  for  43 
points  along  the  .x-direction  and  34  points  in  the 
r-direction.  Since  calculations  are  made  at  incremental 
distances  of  0.1  mm  in  these  two  directions,  the  mapped 
area  is  4.2  mm  by  3.3  mm,  which  means  that  there  are 
1462  calculated  intensity  values  in  the  mapped  area. 

Gridding  software  is  used  to  establish  as  many 
additional  intensity  values  as  desired,  located  between 
the  calculated  data  points,  and  to  create  contour  lines 
connecting  identical  intensity  values.  Software  was  also 
used  to  find  the  maximum  value  of  the  intensity  in  the 
radiated  energy  in  the  water  in  the  mapped  area  which, 
in  the  case  discussed  here,  was  approximately  0.033  of 
the  maximum  value  of  the  intensity  of  the  Gaussian  beam 
incident  at  the  water-solid  interface  of  the  plate. 

Figure  I  shows  the  intensity  distribution  in  the  mapped 
area.  The  conour  lines  are  1  dB  apart,  starting  with  those 
lines  connecting  the  maximum  intensity  points  in  the 
lobes  and  ending  with  those  lines  connecting  intensity 
points  -lOdB  below  that  value.  No  lines  are  drawn  in 
areas  where  the  intensity  is  more  than  10  dB  below  the 
maximum.  All  points  with  intensity  values  above  the 
maximum  in  water  are  connected  with  one  line,  regardless 
of  how  much  higher  the  intensity  value  is  at  those  points, 
compared  with  the  maximum  intensity  in  water.  Clearly 
this  results  in  a  dark  area,  representing  the  plate  where 
the  intensity  is  generally  very  much  higher  than  in  the 
reradiation  pattern  in  the  lobes  in  the  water. 

The  intensity  representation  in  Figure  1  depicts  the 
situation  at  a  given  time  increment,  in  this  case  34  periods 


Figure  1  Calculated  intensity  distribution  of  an  ultrasonic  field  in 
water  created  by  reradiation  of  an  Ai-Lamb  mode  in  an  aluminium 
plate.  Plate  thickness  d  =  0.9  mm  and  parameter /d  =  2.96MHz  mm 
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Figure  2  Stroboscopic  Schlieren  image  of  radiation  lobes  for 
an  aluminium  plate  in  water.  Parameters  are  the  same  as  used  for 
Figure  1 


after  the  start  of  the  generation  of  the  Lamb  wave  in  the 
plate.  Therefore,  the  reradiation  pattern  around  the  end 
of  the  plate  is  well  established,  the  high-intensity  areas 
in  the  water  represent  the  location  of  the  positive  and 
negative  maximum  pressure  amplitudes  in  the  lobes,  i.e. 
the  ^/  locations  of  the  3.29  MHz  wave  in  the  water.  Based 
on  this  particular  appearance  of  the  low  lobes  radiating 
from  the  end  of  the  plate  one  cannot  determine  whether 
the  water  in  the  lobes  are  in  phase  or  180"  out  of  phase 
as  they  leave  the  plate.  This  is  also  the  case  for  a  Schlieren 
image  of  the  ultrasonic  field  at  the  end  of  the  plate,  as 
shown  in  Figure  2. 

This  image  was  obtained  by  creating  a  Schlieren  image 
of  the  ultrasonic  field  where  the  light  source  consisted  of  * 
a  small  HeNe  laser  whose  output  was  modulated  by  a  ,* 
Bragg  cell.  The  time  delay  between  the  generation  of  the  ^ 
leaky  Lamb  mode  and  the  stroboscopic  illumination  of 
the  resulting  radiation  field  could  be  adjusted.  The  image 
presented  here  shows  the  radiation  field  created  by  an 
incident  ‘tone  burst’  consisting  of  four  cycles,  with  the 
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frequency,  mode,  and  plate  dimensions  being  the  same 
as  used  for  the  calculations  (Figure  I).  Also  visible  in 
Figure  2  are  portions  of  the  transmitted  and  reflected 
tone  burst,  consisting  of  long,  straight  wave  fronts.  Their 
interference  with  reradiation  lobes,  can  be  held  to  a 
minimum  when  the  tone  burst  is  short. 

A  comparison  of  Figures  !  and  2  shows  that  it  is 
possible  to  calculate  and  thus  predict  radiation  patterns 
and  lobe  structures  of  plates  immersed  in  a  liquid  and 
vibrating  in  a  given  Lamb  mode.  However,  the  calculation 
can  also  predict  the  relative  phase  of  the  waves  in  the 
radiation  lobes. 

Figure  3  shows  the  amplitude  distribution  in  the  two 
lobes  generated  at  the  narrow  end  of  the  plate  with  only 
the  positive  .x-direction  particle  displacements  considered. 
The  negative-going  half  cycles  of  the  waves  have  been 
suppressed  in  the  graph  and  thus  it  becomes  evident  that 
the  waves  in  the  two  lobes  are  180  out  of  phase.  The 
contour  lines  within  every  section  of  the  waves  indicate 
the  levels  of  amplitude  in  ten  equal  steps  from  the 
maximum  to  zero  amplitude  with  the  lowest  contour 
corresponding  to  the  location  of  the  zero-amplitudes  of 
the  waves.  The  wavefronts  show  an  increasingly  noticeable 
curvature  with  the  circular  pattern  being  centred  at  the 
centre  of  the  narrow  end  of  the  plate. 

Radiation  lobes  from  the  top  and  bottom  surfaces  of 
the  plate  are  visible  in  Figures  I  and  2.  The  points  where 
the  two  lobes  at  the  top  and  the  bottom  of  the  plate 
originate  are  separated  by  one  half-wavelength  of  the 
leaky  A, -Lamb  mode.  Results  obtained  with  the  FDM 
indicate  that  these  lobes  have  intensity  variations  along 


Figure  3  Calculated  wavefront  pattern  for  radiation  front  the 
narrow  plate  end  considering  only  positive  displacement  amplitudes 


Lamb  wave  radiation:  J  -f.  Piet  et  at 


Figure  4  Calculated  wavefront  panern  for  radiation  from  the  top 
and  bottom  surfaces,  considering  only  positive  displacement 
amplitudes 


their  propagation  directions  which  are  cyclic.  The  lobes 
originating  at  the  plate  top  and  bottom  surfaces,  as  shown 
in  Figure  /,  have  high  intensities  at  the  point  of  origin  for 
the  lobes  close  to  the  narrow  end  of  the  plate  while  the 
lobes  originating  farther  away  from  the  end  of  the  plate 
have  a  low  intensity  at  the  point  of  origin.  However,  if 
the  time  increment  between  generation  and  radiation  of 
the  Lamb  mode  is  changed  the  resulting  high  intensity 
areas  in  the  lobes  move  along  the  propagation  direction. 

The  phase  relationship  between  the  wavefronts  within 
the  lobes  originating  on  the  top  and  bottom  surfaces 
as  well  as  the  narrow  end  of  the  plate  can  be  seen  by 
comparing  the  intensity  distribution  in  the  ultrasonic  field 
(Figure  1 )  with  the  amplitude  distribution  ( Figures  3  and 
4)  which  shows  that  the  waves  in  adjacent  lobes  are  180 
out  of  phase.  Figures  3  and  4  reveal  the  phase  information 
which  is  not  available  from  an  observation  of  the 
experimental  results  shown  in  Figure  2,  or  from  the  results 
of  intensity  calculations  shown  in  Figure  I. 
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Standard  dispersion  curves  relating  the  phase  velocity  of  Lamb  waves  to  the  frequency-plate 
thickness  parameter  fd  indicate  that  in  some  cases  symmetrical  and  antisymmetrical 
mode  velocity  curves  cross  each  other.  Viktorov’s  equations  are  used  to  show  that  the  crossing 
points  are  points  where  the  dispersion  curves  are  discontinuous  so  that  no  distinct  Lamb 
mode  exists  for  these  particular  velocity-yi/  combinations.  Procedures  are  given  to  predict  how 
many  such  points  exist  and  where  they  are  located  in  a  given  range  of  fd. 

PACS  numbers:  43.35.Pt 


INTRODUCTION 

Lamb  modes  are  either  symmetric  or  antisymmetric 
and  their  velocity  of  propagation  depends  on  the  frequency 
of  the  mode  /  and  the  thickness  d  of  the  solid  plate  sup¬ 
porting  the  mode  as  well  as  on  the  bulk  wave  velocities  of 
the  material.  If  one  plots  the  phase  velocity  of  the  various 
symmetrical  Lamb  modes  as  a  function  of  the  product  /  -  d 
one  obtains  a  family  of  curves  none  of  which  crosses  an¬ 
other.  The  same  is  true  for  a  plot  of  the  antisymmetrical 
modes.  However,  if  the  two  families  of  curves  are  plotted 
together  there  are  some  instances  where  symmetric  and 
antisymmetric  curves  seem  to  cross,  as  is  evident  from  a 
representative  example  shown  in  Fig.  1  which  applies  to 
aluminum  plates. 

This  implies  that  the  plate  supports  a  symmetrical  and 
an  antisymmetrical  mode  at  the  same  time,  both  modes 
propagating  with  the  same  phase  velocity,  if  the  plate 
thickness  d  and  the  mode  frequency  /  is  such  that  the 
product  fd  is  the  fd  value  of  the  “crossing  point.”  The 
resulting  particle  displacement  within  the  plate  could  then 
be  a  superposition  of  two  different  modes  with  the  vibra¬ 
tional  mode  of  the  plate  not  being  defined  as  a  Lamb  mode. 
If  one  were  to  assume  that  no  pure  Lamb  mode  exists  at 
crossing  points,  the  defining  equations  given  by  Viktorov' 
should  not  be  satisfied  for  those  particular  sets  of  parame¬ 
ters.  Viktorov’s  basic  wave  equations  (1.2)  lead  to  a  4x4 
determinant,  which  in  turn  lead  to  two  characteristic  equa¬ 
tions,  determining  the  eigenvalues  of  the  wave  numbers  of 
the  modes,  given  by  Viktorov  as  Eqs.  (II.4)  and  (II.5). 
One  can  now  use  these  equations  to  calculate  the  value  of 
the  defining  determinant  for  the  values  of  f  d  and  the 
phase  velocity  of  the  Lamb  modes  at  the  crossing  point. 

A  computation  was  performed,  using  double  precision, 
of  the  value  of  the  determinant  for  the  crossing  point  lo¬ 
cated  near  the  value  of  f  d=S.5  MHz  mm  and  the  phase 
velocity  of  around  6.9396  km/s,  shown  in  Fig.  1.  These 
calculations  for  the  value  of  the  determinant  were  per¬ 
formed  by  changing  the  values  of  f  d  in  increments  of 
10“*  MHz  mm  and  c  in  increments  of  0.1  mm/s.  The 
result  is  that  for  the  velocity  range  c=6.9396744  to 
c= 6.9396764  km/s  and  the  f  d  range  from  8.45435146  to 
8.45435157  MHz  mm  the  determinant  is  nonzero.  This  re¬ 
sult  would  indicate  that  a  Lamb  mode  does  not  exist  at 


points  where  the  velocity  dispersion  curves  for  the  sym¬ 
metrical  and  the  antisymmetrical  modes  cross. 

The  calculations  leading  to  this  result  are  rather  in¬ 
volved  and  would  be  very  taxing  if  all  the  crossing  points 
had  to  be  found  with  no  prior  knowledge  of  the  approxi¬ 
mate  location  of  these  points.  Fortunately,  a  much  shorter 
method  is  available  to  determine  the  location  of  those 
points  on  the  families  of  velocity  dispersion  curves  where 
no  solution  to  the  defining  Lamb  equations  exist.  This 
method  will  be  discussed  below. 


I.  THEORETICAL  CONSIDERATIONS 


Noting  that  the  phase  velocities  of  Lamb  modes  at 
critical  crossing  points  are  greater  than  the  bulk  longitudi¬ 
nal  wave  velocity,  one  can  use  the  defining  Lamb  mode 
equations  for  the  symmetrical  and  antisymmetrical  modes, 
respectively,  as  given  by  Viktorov* 

tan  a/tan  P+  B=0, 

(1) 

tan  /3/tan  a+B=0, 

(2) 

where 

a=(npfd)/Cp 

(3) 

P={irqfd)/c„ 

(4) 

p=  Vl-(c/c)^ 

(5) 

q=  ^/(c,/c,)^-(c/c)^ 

(6) 

4(cA)V<7 

^  [2(c/c)‘-l] 

(7) 

with  c,  the  shear  wave  velocity  of  the  bulk  solid,  C/  the 
longitudinal  wave  velocity,  and  c  the  phase  velocity  of  a 
Lamb  mode  defined  by  Eqs.  (1)  and  (2). 

As  written,  Eqs.  (1)  and  (2)  apply  as  long  as  c>c;. 
This  is  the  region  of  interest  for  the  present  purpose.  One 
can  write  general  expressions  for  the  defining  equations 

5(c,/rf)=tana/tan/8-f  5,  (8) 

Aic,fd)  =tiin  P/tan  n  + B.  (9) 
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FIG.  1.  Section  of  Lamb  mode  velocity  dispersion  curves  for  aluminum 
showing  the  phase  velocity  c  and  thickness-frequency  parameter  fd  for  the 
first  four  mode  crossing  points. 

Solutions  of  Viktorov’s  equations  are  determined  by  the 
roots  of  S{c,fd)  and  A(c,fd),  i.e.,  the  values  of  c  and  an 
appropriate  fd  are  defined  for  an  S  mode  and  an  A  mode 
when  S{c,fd)  or  A{c,fd),  respectively,  are  zero. 

II.  CRITICAL  POINTS  OF  VIKTOROV’S  EQUATIONS 

The  variables  a  and  0  are  the  functions  of  the  param¬ 
eters  c  and  fd.  Thus,  for  any  given  phase  velocity  c,  tan  a 
and  tan  P  are  periodic  functions  of  fd,  where  the  periodic¬ 
ities  are  cfp  and  cfq,  respectively.  The  arguments  of  both 
tangent  functions  are  zero  for  fd=Q,  regardless  of  the  value 
of  c.  With  increasing  values  of  fd,  the  values  of  S(c,fd) 
and  Aic,fd)  usually  remain  nonzero;  however,  in  order  to 
satisfy  Viktorov’s  equations  they  must  be  zero  which  oc¬ 
curs  only  when  either  tana/tan/3=— 5  or  tan /3/tan  a 
=  —B  according  to  Eqs.  (8)  and  (9). 

These  conditions  on  Eqs.  (8)  and  (9)  must  be  satisfied 
simultaneously  at  mode  crossing  points  in  the  c-fd  plane  if 
an  A  mode  as  well  as  an  5  mode  are  to  exist  at  such  points. 
Investigating  the  behavior  of  the  two  tangent  functions  one 
finds  that  two  types  of  critical  points  may  exist  where  there 
are  no  solutions  to  Viktorov’s  equations.  As  will  be  shown 
below,  these  critical  points  in  the  c-fd  plane  are  identical 
to  the  mode  crossing  points;  one  of  these  points  was  iden¬ 
tified  in  the  introductory  discussion  of  the  determinant  so¬ 
lution. 

Type  I  critical  points  exist  when  the  respective  values 
of  tan  a  and  tan  /3  both  reach  infinity  at  the  same  value  of 
fd  which  occurs  when 


a=  (.irpfd)/c,=  im-\-{)7r. 

(10) 

P={'trqfd)/c,=  (n-\-{)Tr, 

(11) 

for  certain  values  of  the  integers  m  and  n  and  for  certain 
values  of  p  and  q  which  in  turn  are  functions  of  the  phase 

velocity. 

Type  II  critical  points  exist 

when  tan  a  and  tan  P 

reach  zero  at  the  same  value  of  fd  that  requires  that  the 

arguments  are 

a={Trpfd)/c,=mir, 

(12) 

P=(Trqfd)/c,—niT. 

(13) 

The  procedure  to  determine  at  what  values  of  c,  fd, 
and  m  and  n  type  I  and  type  II  critical  points  occur  for  a 
set  of  Lamb  mode  velocity  dispersion  curves  consists  of 
combining  Eqs.  (10)  and  (11)  for  the  type  1  points,  and.' 
Eqs.  (12)  and  (13)  for  type  II  critical  points.  •. 

As  an  example  consider  the  procedure  for  type  I  crit¬ 
ical  point  evaluation  for  an  aluminum  plate  where  the  lon¬ 
gitudinal  and  shear  bulk  wave  velocities  are,  respectively, 
6.42  and  3.04  km/s.  These  are  the  same  values  that  were 
used  in  the  determinant  calculation  mentioned  in  the  In¬ 
troduction  above. 

Treating  m  and  n  as  integers  whose  magnitude  is  not 
specifically  assigned,  Eqs.  ( 10)  and  ( 1 1 )  are  two  equations 
with  two  unknowns,  i.e.,  a  critical  value  of  fd  and  a  critical 
value  of  c  in  and  q,  respectively.  Not  specifying  for  the 
moment  the  value  of  fd  but  knowing  that  fd  will  have  to  be 
the  same  in  a  and  /3  where  a  critical  point  is  located,  one 
can  substitute  the  expression  for  fd  from  Eq.  ( 10)  into  Eq. 
(11)  and  solve  for  c  in  terms  of  m  and  n.  This  yields  the 
following  expression  for  the  critical  phase  velocity,  cy. 

I 

Um-|-i)^(c/c,)^-(n-f  j)^ 

For  the  critical  phase  velocity  to  be  real,  nonzero,  and 
greater  than  the  longitudinal  wave  velocity,  the  conditions 
for  m  and  n  are 

m>n,  (m-(-5)/(n-f5)  >c//c, =2.12  for  Al. 

For  these  conditions  and  the  condition  that  tan  a  and  tan 
P  go  to  infinity  at  the  same  fd,  one  finds  only  three  (m,n) 
combinations,  (1,0),  (2,0),  and  (3,0)  for  the  range  of  the 
parameter  fd  given  by  0<yif  <  12. 

Substituting  these  three  sets  of  values  for  (m,n)  into 
Eq.  (14)  yields  three  critical  phase  velocities.  The  three 
corresponding  critical  fd  values  are  found  by  substituting 
the  three  sets  of  and  (m,n)  into  either  Eq.  (10)  or  Eq. 
(11)  to  solve  for  the  value  of  fd.  The  three  critical  points 
for  aluminum  are  thus  found  at 

c=8.5221  km/s,  /rf=4.88  MHz  mm,  (m,n)  =  (l,0), 

c=6.9396  km/s,  /d=8.45  MHz  mm,  (w2,n)  =  (2,0), 

c= 6.6647  km/s,  /</=  11.96  MHz  mm,  (m,n)=3,0). 

Type  II  critical  points  are  determined  in  a  similar  fash¬ 
ion  by  using  Eqs.  (12)  and  (13).  Following  essentially  the 
same  procedure  as  outlined  above,  one  finds  the  possible 
critical  velocities  to  be  given  by 

I 

Making  the  appropriate  substitutions  one  finds  only  one 
(m,n)  set  in  the  range  0  <yif  <  12.  This  single  type  II  crit¬ 
ical  point  has  the  values 

c= 8.5221  km/s,  /d= 9.76  MHz  mm,  (m,«)  =  (3,l). 

Comparing  these  results  with  the  mode  velocity  curves 
of  Fig.  1  shows  that  the  four  calculated  critical  points  are 
the  only  crossover  points  of  A  and  5  modes  in  the  range 
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0<fd<\2.  Inasmuch  as  the  parameters  of  the  critical 
points  do  not  satisfy  Viktorov’s  defining  equations,  the 
'  crossover  points  are  not  valid  solutions  either  and  neither 
*Tan  antisymmetrical  nor  a  symmetrical  mode  exists  at  these 
^points. 

III.  CONCLUSION 

The  Lamb  mode  velocity  dispersion  curves  generally 
exhibit  a  number  of  points  where  an  A  mode  and  an  S 
mode  cross  which  seems  to  imply  that  for  that  particular 
set  of  parameters  a  solid  plate  should  be  able  to  support  a 
symmetrical  as  well  as  an  antisymmetrical  mode.  Examin¬ 
ing  the  simplest  form  of  Viktorov’s  equation  shows  that  the 
parameters  defining  these  crossing  points  do  not  constitute 


a  solution  to  the  equations,  thus  these  unique  points  are  not 
parts  of  either  the  symmetrical  or  antisymmetrical  velocity 
dispersion  curves.  The  relatively  simple  procedure  de¬ 
scribed  here  can  be  used  to  find  the  number  of  crossing 
points  of  Lamb  mode  velocity  dispersion  curves  as  well  as 
their  velocity  and  fd  parameters  of  these  points. 
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